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ABSTRACT 

We investigate 3D atmosphere dynamics for tidally locked terrestrial planets with 
an Earth-like atmosphere and irradiation for different rotation periods {Prot = 1 — 
100 days) and planet sizes {Rp = 1 — ^Rparth) with unprecedented fine detail. We 
could precisely identify three climate state transition regions that are associated with 
phase transitions in standing tropical and extra tropical Rossby waves. 

We confirm that the climate on fast rotating planets may assume multiple states 
{Prot ^ 12 days for Rp = 2REarth)- Our study is, however, the first to identify the type 
of planetary wave associated with different climate states: The first state is dominated 
by standing tropical Rossby waves with fast equatorial superrotation. The second state 
is dominated by standing extra tropical Rossby waves with high latitude westerly jets 
with slower wind speeds. For very fast rotations {Prot ^ 5 days for Rp = 2REarth), we 
find another climate state transition, where the standing tropical and extra tropical 
Rossby wave can both fit on the planet. Thus, a third state with a mixture of the two 
planetary waves becomes possible that exhibits three jets. 

Different climate states may be observable, because the upper atmosphere’s hot 
spot is eastward shifted with respect to the substellar point in the first state, westward 
shifted in the second state and the third state shows a longitudinal ’smearing’ of the 
spot across the substellar point. 

We show, furthermore, that the largest fast rotating planet in our study exhibits 
atmosphere features known from hot Jupiters like fast equatorial superrotation and a 
temperature chevron in the upper atmosphere. 

Key words: planets and satellites: atmospheres -planets and satellites: terrestrial 
planets - methods: numerical. 


1 INTRODUCTION 

Transiting habitable Super-Earths around small, cool 
M dwarf stars are easier to detect than other habitable 
planet-star configurations of the same stellar brightness 
due to their favourable star-p lanet contrast ratio (e.g. 
iMadhusudhan fc Redfieldl ll2015l 'll. They are thus the prime 
targets for the characterization of the atmospheres of terres¬ 
trial habitable exoplanets in the near future. These planets 
allow, furthermore, to test planet formation, evolution and 
planet atmosphere properties under conditions, for which 
there is no counterpart in the Solar System. Several up - 
coming spa ce missions, such a s TESS (iRicker et al.l 120141 , 
CHE OPS llBroeg et al.l l2013ll a nd PLATO llRauer et all 
2 OI 4 II. and g r ound based surveys llSnellen et al.l (l2012tl and 
Gillon et all (120131 '). among others) are expected to yield 
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candidates for the characterization of Super-Earths planets 
around bright nearby stars. 

At least the planets at the inner edge of the habit¬ 
able zone of M dwarfs are so close to their star that tidal 
interactions can force the planet’s rotation in short astro¬ 
nomical time scales to be syn c hronous with the orbital p e- 
riod le.g. iKasting et aP (ll993l 'l. lMurrav fc DermottI ll200dl 'l. 
They should thus continuously face their star with the same 
side; resulting in an eternal day and night side. Such a 
configuration is definitely not Earth-like and challenges our 
concept of ’habitability’. Consequently, there are still many 
uncertainties about the exact shape of planet atmosphere 
dynamics that such alien planets can assume. Parameter 
studies are one tool to investigate different climate dynamic 
regimes on tidally locked planets, but are very costly in com¬ 
putation time. 

Ijoshi et al.l (If 9971 ') were the first to investigate 3D large 
scale planet atmosphere dynamics on tidally locked habit- 
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able planets for selected parameters. They found that heat 
can be efficiently transported from the day side towards 
the night side for sufficiently dense atmospheres (surface 
pressure ps ^ 100 mbar). They focused, however, on small 
Earth-size planets and did not change the rotation period. 

Several more parameter studies have been performed in 
the last years that investigated various aspects of 3D planet 
atmosphere dynamics on tida lly locked habitabl e planets 
with an Earth-like atmosphere: lEdson et al.l |20l3) revealed 
an abrupt phase change in atmosphere dynamics on an 
Earth-size planet fo r rotation periods between 3 — 5 days. 
iHeng fc Vogd (|201ll l investigated surface temperatures and 
wind structures for different surface friction, radiative time 
scale assumptions, surface gravity, surface pressures and ro¬ 
tation periods on a Super-Earth planet {Rp = 1.45REarth)- 
They showed that different surface friction assumptions lead 
to very different su r face w ind and temperature structures. 
iMerlis fc Schneider! (120101 1 investigated the general climate 
state for a tidally locked Earth-size ocean planet for very 
fast {Prot = 1 day) and very slow rotation {Prot = 365 days) 
and found that the atmosphere dynamics of the fast rotator 
is dominated by westerly jets at high latitudes, whereas the 
dynamics on the slow rotator is dominated by direct circula¬ 
tion fr om the substellar point to the night side. lCarone et aP 
||2014| investigated the dynamics of a Super-Earth planet 
[Rp = lA^REarth) for two rotation periods Prot = 10 and 
Prot = 36.5 days and showed, among other things, the ex¬ 
istence of an embedded reverse circulation cell inside in a 
larger direction circulation cell for the faster rotation. 

These investigations used models of different complex¬ 
ities and water content, focused on different aspect of at¬ 
mosphere dynamics, and - what is more important - cov¬ 
ered only a small part of the relevant planet rotation and 
planet size parameter space for tidally locked habitable plan¬ 
ets around M dwarfs. The current state of knowledge makes 
it thus difficult to draw conclusions about the evolution of 
climate dynamics with increasing rotation. Even worse, cli¬ 
mate state changes can be overlooked in a ’coarse’ study, if 
they appear unexpectedly in a small part of the parameter 
space. It is, furthermore, difficult to identify the underly¬ 
ing mechanisms for climate state transitions as the usage of 
different planet parameters and model prescriptions act as 
confounding factors. 

The purpose of this paper is, thus, to remedy this short¬ 
coming: We will map the climate dynamic regime for terres¬ 
trial habitable planets around M dwarfs with hitherto un¬ 
precedented thoroughness, covering Prot = 1 — 100 days and 
Rp = 1 — 2REarth- We will demonstrate that the pertur¬ 
bation method and the meridional Rossby wave numbers 
constitute very powerful diagnostic tools for the discussion 
of the simulation results. The two methods facilitate the 
precise identification of climate phase changes and multi¬ 
ple climate solutions for the same set of parameters. They 
allow us, furthermore, to identify the specific Rossby wave 
responsible for a given climate state. We can even coherently 
compare climate dynamics arising from different models and 
explain the likely origin of discrepancies. 

In the following, we will first briefly describe the model 
and numerical adjustments. We will then introduce useful 
quantities like tropical and extra tropical Rossby radius of 
deformation, dynamical and radiative time scales, and the 
perturbation method. We will use these quantities to identify 


climate regime transitions due to different standing Rossby 
waves in the investigated parameter space. 

With an improved understanding of the Rossby wave 
climate regimes, we will then discuss the evolution of differ¬ 
ent atmosphere features like zonal wind jets, cyclonic vor¬ 
tices, the state of circulation cells, day side and night side 
surface temperatures, and vertical temperature structure. 
We will identify characteristic features that are associated 
with different standing Rossby wave regimes: the location of 
the upper atmosphere hot spot, the number and strength of 
zonal wind jets, and the state of circulation cells. 

We will then briefly discuss the potential to observation- 
ally discriminate between different climate dynamic states, 
provide a summary of our results and finish with a conclu¬ 
sion and outlook. 


2 THE MODEL 

Our model is described in details in ICarone et al.l ll20l4l . 
It is suitable for the 3D modelling of tidally locked terres¬ 
trial planets with zero obliquity. It can be used for different 
atmosphere compositions, planetary and stellar parameters 
(stellar effective temperature and stellar radius). The model 
was developed to be as versatile and computationally effi¬ 
cient as possible, using simple principles and the atmosphere 
properties of Solar System planets as a guide line. We will 
give here a brief prescription and explain the specific set-up 
used for this particular study. 

2.1 The dynamical core: MITgcm 

We employ the dynamical core of the Massachusetts Insti¬ 
tute of Technol(^y global circulation m odel (MITgcm) de¬ 
veloped at MldJj (lAdcroft et al.l l2004f l. It uses the finite- 
volume method to solve the primitive hydrostatical equa¬ 
tions (HPE) can be written as the horizontal momentum 
equation 


Dv -* 

— + fk XV = -Vp® - 1 - Tv, 

(1) 

vertical stratification equation 

d-l- 1 

dp p’ 

(2) 

continuity of mass 

^ 9uj 

Vp ■ V + — = 0, 
op 

(3) 

equation of state for an ideal gas 

p = p{R/p) T, 

(4) 

and thermal forcing equation 

DO e 

'm ^ 

(5) 

where the total derivative is 

D d ^ „ d 

(6) 


^ http://mitgcm.org 
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Here, v is the horizontal velocity and oj = Dp/Dt is the ver¬ 
tical velocity component in an isobaric coordinate system 
and, p and t are pressure and time, respectively. ^ = gz 
is the geopotential, with 2 being the height of the pressure 
surface (also called geopotential height), and g the surface 
gravity. / = sin v is the Coriolis parameter, where is 
the planet’s angular velocity and v the latitude at a given 
location, k is the local vertical unit vector, T, p and Cp are 
the temperature, density, and specihc heat at constant pres¬ 
sure, respectively. 9 = T (p/pa)'^ is the potential tempera¬ 
ture with K being the ratio of the specihc gas constant R/p 
to Cp. Here, R = 8.314 JK~^mol“^ is the ideal gas constant, 
p is the molecular mean mass of the atmosphere, and Pa is 
the mean constant surface pressure dehned at the surface. 
The pressure at surface level piQmbar) varies horizontally 
around Pa in the model. Re,v are the horizontal momentum 
and thermal forcing terms, respectively. Note that is a 
vector and Te is a scalar. 


2.2 Temperature forcing, surface friction and 
sponge layer 


The temperature forcing Re is realized by using Newtonian 
cooling that drives the atmosphere temperature T with a 
relaxation timescale T^ad towards a prescribed equilibrium 
temperature T^q, which writes in the easier recognizable 
temperature form as: 




( 7 ) 


Trad is Calculat ed us ing basic estimates (e.g. 
IShowman fc Guillotl (l2002l 'l'). see also Section [3|. The 
equilibrium t e mper ature T^q is described in details in 
ICarone et al.l (|2014h . where we have developed Teq for 
the day and night side. The illuminated dayside is driven 
towards radiative-convective temperature profiles suitable 
for a greenhouse atmosphere. The nightside temperature 
relaxes towards the Clausius-Clayperon relation of the main 
constituent of the atmosphere. Furthermore, we capped 
the top of the atmosphere with a vertically isothermal 
stratosphere. 

The horizontal momentum forcing = Rfric + Raponge 

consists of a surface friction and a sponge layer term: 


R fric 


fric 


( 8 ) 


and 


R SI 


= —kRV, 


( 9 ) 


where we use the Rayleigh friction approximation in both 
cases. The first term acts at the surface in the planetary 
boundary layer (PEL) between pa and 0.7 x pa- The sur¬ 
face friction time scale Tsjric is set to 1 day and linearly 
decreases with pressure until the top of the PEL is reached. 
The choices in boundary layer extent and surface fric tion 
time scale are taken over from iHeld fc Suared l|l994|) for 
the E arth. However, we already discussed in ICarone et ^ 
(l2014h that the prescription of the boundary layer varies 
greatly between models, even for the best studied other ter¬ 
restrial planet: Mars. We will investigate changes in bound¬ 
ary layer extent and their effect on climate regime for tidally 
locked terrestrial paper in an upcoming paper. For now, we 


Table 1. Sponge layer friction timescales used in our model 


k 

•^max 

Prot 

[day-1] 

[days] 

1/80 

4-100 

1/40 

2-4 

1/20 

^ 2 


keep the Earth values from iHeld fc Snared lll994h as de¬ 
scribed here. 

The sponge layer term stabilizes the upper atmosphere 
boundary layer numerically against non-physical wave reflec¬ 
tion. It may also be justihed physically as a crude represen¬ 
tation of gravity wave breaking at low atmosphere pressure. 
The sponge layer is implemented for p ^ 100 mbar, where kn 
increases from zero to kmax at t he upper atmosphere b ound ¬ 
ary using the pres cription of IPolvani fc Kushnerl (l2002li 
(ICarone et al.l 120141 '). In the nominal set-up kmax was set 
to 1/80 day”*^. However, we found that the sponge layer ef- 
hciency needs to be increased in the fast rotation regime to 
ensure numerical stability. It will be shown later that this 
regime exhibits very fast equatorial superrotating jets with 
wind speeds larger than 100 m/s at the top of the atmo¬ 
sphere that suppresses upwelling over the substellar point. 
The dynamical interplay between these two dynamical fea¬ 
tures triggers more vigorous upper atmosphere dynamics 
that requires a stronger sponge layer to prevent un-physical 
wave reflection at the upper boundary of the model. The 
adjustment of sponge layer friction may tentatively be jus¬ 
tified physically by assuming a more efficient wind braking 
mechanism for very fast wind speeds. See Table [T] for the 
adopted kmax dependent on rotation period. Conventionally, 
the sponge layer is omitted when discussing results because 
the layer is assumed to be not physically ’real’. 

It should be noted, though, that Rt and R/ric are ideal¬ 
ized representations of external heating and boundary layer 
physics, respectively. 


2.3 Numerical setup 


We use the same set-up as described in ICarone et al.l ll2014h . 
We use 20 vertical levels with equally spaced 50 mbar levels 
and surface pressure Pa = 1000 mbar, to resolve the tro¬ 
posphere, that is, the w eather layer of a nd E arth-like at¬ 
mosphere as outlined in Held_J^_Suarej (119941 1. Indeed, it 


was confirmed by IWilliamson et ahT i 199^1 that this vertl 


cal resolution is sufficient to resolve the general structure 
of the Earths troposphere. This vertical resolution is conse- 


Earth-like planets (iHen 

e k, Vofft 20 111: Edson et al. 

201 1 |; 

iMerlis & Schneideij|201C 

: Joshi et al. 1997: IJoshil 2003 

). We 


use the C32 c ubed-sphere grid for solving in the horizon¬ 
tal dimension ([Marshall et al.ll2004l l. which corresponds to 
a global resolution in longitude-latitude of 128 x 64 or 
2.8° X 2.8°. The substellar point is placed at 0 degrees lon¬ 
gitude and latitude. _ _ 

As established in ICarone et al.l (l2014h . we evolve the 
atmosphere from initial conditions of universally constant 
temperature T = 264 K for tinu — 400 Earth days. Data 
generated before tinu are discarded; the simulation is run 
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subsequently for trun = 1000 Earth days and averaged over 
this time period. For Rp ^ 1.25REar th, the nominal tim e 
step At = 450 s is applied, as used in ICarone et al.l l|2014h . 
For the Earth-size planet simulations, At = 300 s is required 
to maintain the Courant - Friedrichs -Lewy (CFL) condition 
for the numerical stability of the finite volume method used 
in MITgcm. A reduction in planet size, while maintaining 
the C32 cube sphere grid, reduces the horizontal length scale 
and thus necessitates the reduction in time step. 


2.4 Investigated parameter space: Tidally locked 
habitable terrestrial planets aronnd M dwarfs 

We aim to investigate climate dynamics on tidally locked 
habitable planets around M dwarf stars more thoroughly 
than previous studies with respect to rotation period and 
planet size. At the same time, we strive to minimize con¬ 
founding factors by keeping, e.g., thermal forcing constant 
across the investigated parameter space. The combination of 
a very thorough parameter sweep and the minimization of 
confounding factors allows to identify more climate dynam¬ 
ics phase changes and to better understand the underlying 
mechanisms behind them. It will be confirmed that in par¬ 
ticular the rotation period is of great importance for the 
evolution of climate dynamics. 

We keep for this study an Earth-like atmosphere com¬ 
position: A nitrogen dominated greenhouse atmosphere with 
surface pressure pa = 1000 mbar. See Table[2]for the relevant 
parameters 0. We assume, furthermore, that each planet re¬ 
ceives the same amount of net stellar energy than the Earth. 
Such a planet lies in any case in t he ha b itable zone as in¬ 
vestig a ted by, e.g., iKasting et al.l lll993 ): iKopparapu et ^ 
(l201^ : lYang et all ( 2014 1: Zsom et al.l I 2013f b 

We introduce the net incident stellar flux, because not 
all of the incident stellar energy is used to heat a planet’s 
atmosphere; part of it is reflected back, e.g., by clouds and 
ice. In the grey two-stream greenhouse model used in this 
study to calculate radiative-convective equilibrium temper¬ 
atures Teq, we encapsulate the amount of reflected stellar 
energy in the planetary albedo a. The net incoming stellar 
flux is thus 


fne£ — fo(l n;). 


( 10 ) 


where the incident stellar flux that the planet receives 
at semi major axis a is: 


lo 


La 

47ra^ ’ 


( 11 ) 


where L* is the stellar luminosity and a can be calculated 
from the orbital period and stellar mass using Kepler’s third 
law. For the Earth, the following values can be assumed Jo = 
1368 W/m~^ and a = 0.3, which yields Rat = 958 W/m“^ 
dCarone et al.l[20l3 l. 

This study is highly relevant in mapping the climate dy¬ 
namics landscape in preparation for other GCMs that self- 
consistently calculate a by taking into account cloud and 
surface ice formation, but are too computationally costly 
to perform a detailed parameter sweep in orbital period. 



Figure 1. Orbital period for tidally locked planets around M 
main sequence stars receiving Ret = 958 W/m“^ for planetary 
albedo a = 0,0.3, 0.5, 0.7, 0.9. 


The usage of a fixed Ret in the parameter study allows us 
- in contrast to other GCMs - to avoid changing a, which 
would change thermal forcing. The elimination of changes in 
the thermal forcing as a confounding factor brings the de¬ 
pendence of climate dynamic regimes on rotation period to 
the fore. We can thus identify particularly interesting rota¬ 
tion periods that are ’turning points’ for climate dynamics. 
We will explore in the following, for which orbital period 
range around M dwarfs it is indeed reasonable to assume 
Ret = 958 W/m“^. 

To derive luminosity and mass of M dwarfs, we use Yale- 
Yonsei isochrones for the low stellar mass regime with ap¬ 
proximately solar values: an stellar age of 4 Gyrs, hydro¬ 
gen mass fraction X = 0.70952, Z = 0.01631, metallicity 
\Fe/H] = 0.0 and mix ing length scaled o ver pressure scale 
height umlt = 1.875 (ISoada et al.ll20l3l . The Yale-Yonsei 
isochrones cover M dwarf spectral cl asses M0V-M6V with 
masse s between M* = 0.1 — 0.6M.c;„r,. lKaltenegger fc TraubI 
ll2009ll list smaller stars in their study but note that their 
masses are not well constrained. 

We can then calculate with Equations and (HD), 
the orbital and thus rotation periods of tidally locked plan¬ 
ets around M dwarfs that receive Ret = 958 W/m~^, as¬ 
suming 0 = 0, 0.3, 0.5, 0.7, 0.9. We find that Earth-like ther¬ 
mal forcing is indeed possible in the orbital period range 
Prot = 1 — 100 days for M dwarf starfl if we a llow e xtreme 
values of a ('Figure [Til. iKaltenegger fc TraubI (l2009ll derive 
Prot = 2 — 65.5 days by extending their M dwarf range to 
MOV stars, which is not so different from our rotation period 
range, but where the authors assumed a constant a = 0.3. 
While it is beyond the scope of this work to decide which 
planetary albedos are reasonable for tidally locked terres¬ 
trial p lanets around M dwarfs, we note that iLeconte et al.l 
(l2013ll consider a = 0.1 — 0.5, which corresponds to a rel¬ 
evant orbital_geriod_range between Porb = 6 — 80 days. 
Also. lYang et al.l (l2013h derive a planetary albedo as high as 
a = 0.5 for tidally locked planets around M dwarfs due to 


^ Note that we use R p = lAbREar t h inst ead of Rp = 1-^REarth 
to be consistent with ICarone et alj ll2014ll . 


^ Older/younger M main sequence stars are more/less luminous 
than the stars investigated here. The corresponding planet orbits 
would thus shift towards/away from the star, respectively. 
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Table 2. Planetary and atmospheric parameters used in this 
study for all simulations 


Net incoming stellar flux Inet 

958 W/m2 

Obliquity 

0 ° 

Cp 

1.04 J/gK 

Vs 

1000 mbar 

Main constituent 

N 2 

Molecular mass fi 

28 g/mol 

Atmospheric emissivity e 

0.8 

Surface optical depth Ts 

0.62 

Adiabatic index 7 

7/5 

Surface friction timescale Tgjric 

1 day 

Rotation period Prot 

1-100 days 


Table 3. Planetary size, surface gravity, and mass assuming 
Earth-like bulk density 


Rp 

[PEarth] 

g 

[m/s^] 

Mp 

[Mparth] 

1 

9.8 

1.0 

1.25 

12.5 

2.0 

1.45 

14.3 

3.1 

1.75 

17.3 

5.4 

2 

19.6 

8.0 


stabilizing cloud feedback above the substellar point, which 
extends the inner edge of the habitable zone towards the 
star. 

Fast rotation periods with Prot ^ 6 days allow us as a 
’bonus’ to study dynamical features already known from hot 
Jupiters. In particular, the largest planet size investigated 
in this study {Rp = 2REarth) will be of interest to bridge 
the gap between terrestrial planets and hot Neptunes. In 
other words, this study already lays the foundation for a 
coherent atmosphere dynamics study across the transition 
region between rocky and gas giant planets. 

The selected rotation perio d range Pro t = 1 — 

100 days was also i nvesti gated by lEdson et al.l ll201ll l and 
iMerlis fc Schneider! ll2010l l. The authors focused, however, 
on an E arth-size planet, w hereas we also investigate larger 
planets. lEdson et al.l ((2oT3) scanned the whole rotation pe¬ 
riod region Prot = 1 — 100 days sparsely and focused on a 
specific rotation period region between Prot = 3 — 5 days, 
which the authors associated with climate bifurcation. That 
is, the authors inferred the existence of two different at¬ 
mosphere dynamic states for the exac t same set of plan¬ 
etary parameters and rotation periods. iMerlis fc Schneider! 
ll2010li investigated two rotation periods, Prot = 1 days and 
Prot = 365 days. The latter case is dynamically very si milar 
to the Prot = 100 days simulation of lEdson et al.l (1201 il l, the 
dynamics of which is dominated by direct circulation flow. 
We will show later that direct circulation is expected to dom¬ 
inate climate dynamics for s l ow rot ators {Prot >>45 days). 

Recently. iLeconte et al.l ll2015l ') questioned the common 
assumption that planets in the habitable zone of an M dwarf 
are tidally locked. However, the climate of tidally locked 
planets is very different from asynchronous rotating planets 
like the Earth. It gives rise to new interesting climate dy¬ 
namics states, as we will show in this study. Eurthermore, 
planets at the inner edge of the habitable zone are still prob¬ 


ably tidal ly locked - even in t he light of the revised tidal 
theory of ILeconte et al.l (l2015l l. Incidentally, these planets 
should be the first that become accessible to atmosphere 
characterization. 

We investigate in this study planets of terrestrial com¬ 
position with planet sizes between Rp = 1 — 2REarth- 
This planet size range covers the Super-Earth planets, 
for which rocky composit ion can be safely assumed (e.g. 
iMadhusudhan et al.l 1120121! '). In addition, we assume Earth¬ 
like bulk density, pEarth = 5.5 g/cm®, for easier comparison 
between radiative and dynamical time scales and thus easier 
analysis of atmosphere dynamics evolution with planet size. 
Surface gravity can then be calculated via 


g = 


GMp 


( 12 ) 


where G is the gravitational constant and planetary mass 
Mp can be derived by: 

Mp = pEarth—ltF^P, (13) 

which yields planetary masses between 1-8 Earth masses 
(Table El). 

We are aware that we probably underestimate the mass 
and thus the surface gravit y of large roc ky terrestrial Su¬ 
per Earth planets (see e.g. iHatzed ll2013l! !. We emphasize 
again, however, that we strive in this work to identify cli¬ 
mate dynamics tendencies with rotation period and planet 
size. Thus, it is justified to minimize other confounding fac¬ 
tors like differences in mean density between large and small 
terrestrial planets. Furthermore, there are still many un- 
certainties in the composit i on of Super-Earth plane ts (e.g. 
IMadhusudhan et al.l (120121! . [Rogers fc Seagetl (120101! !. 


3 ATMOSPHERE DYNAMICS TIME AND 
LENGTH SCALES 


For the discussion of climate state evolution with respect 
to rotation period and evolution, we now introduce basic 
atmosphere dynamics time and length scales. The most 
important time scales in the investigation of irradiated 
planet atmospheres are the dynamical Tdyn and the radia- 
tive time scale Trad- T he former is for tidally locked planets 
(IShowman et al.l|201ll! 

Tdyr^ OC (14) 

where U is the mean horizontal velocity or the horizontal 
velocity scale and the length scale L is set to Rp, because 
circumglobal flow driven by the longitudinal heat gradient 
between the substellar point and th e night side is assum ed. 
The radiative time scale, as used in ICarone et al.l (I2ni4h . is 
defined as 


CpPs 

AgaPr 


(15) 


where Ta is the local surface temperatu re and a is the S tefan - 
Boltzmann constant. We discussed in ICarone et al.l (l2014l ! 
the usefulness and limitation of this simple estimate and 
concluded that it is suitable for the modelling of the shallow 
atmosphere of terrestrial planets. 
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For the investigation of atmosphere dynamics with re¬ 
spect to planet rotation period, the most important param¬ 
eter to monitor is the Rossby radius of deformation. Rossby 
waves are pla netary waves that conserve absolute vorticity 
llHoltonlll992fl V ~ C + where (^ = fc-Vxu = is the (verti¬ 
cal component of the) relative vorticity and / = 2fl sin v is 
the Coriolis parameter. Rossby waves owe their existence to 
variations of the Coriolis parameter with latitude and typ¬ 
ically propagate in westward direction. The Rossby radius 
of deformation represents the maximum amplitude of the 
planetary waves that can be excited at a given latitude. 

The extr a tropical Rossby radius is dehned as (e.g. 
iHoltorJ lll992h l 


Lr = 


NH 


(16) 


where N = \J —pg^ fQ x AO/dp is the buoyancy frequency in 
pressure coordinates that encapsulates the static stability 
of the atmosphere and H = RT/pg is the scale height. We 
hnd the best agreement with changes in planet atmosphere 
dynamics, if we evaluate Lr for v = 30°. 

The tropical Rossby radius is dehned as 




NH 


(17) 


where /3 = 2Q./Rp i n the equatorial beta plane approxi¬ 
mation (iHoltonI I 1992 I I. We compare atmosphere dynamics 
for different ratios of the Rossby radius of deformation over 
the planet’s radius Rp. It should be noted that we use 
here a different dehni t ion fo r Rossby radius of deformation 
th an in Carone et al.| |2014[l. where we used the dehnition 
of IShowman fc Polvani ( 20111 1 for a shallow water model. It 
will be shown in this study that the dehnitions presented 
above, that are also used bv lEdson et al.l (l201ll l. are more 
appropriate for a fully 3D global circulation model. 

We will show that the meridional Rossby wave numbers, 
that is, the ratios of the Rossby radii of deformation over the 
planet’s radius Rp, Lr/Rp and Xr/Rp, are extremely im¬ 
portant diagnostic tools in identifying possible climate state 
transitions with respect to rotation period and planet size. 
We will, thus, use these properties extensively to classify and 
discuss climate states. 


4 THE PERTURBATION METHOD FOR 
ROSSBY WAVE ANALYSIS 

For a better qualitati ve analysis of atm osphe ric waves, we 
follow the example of lEdson et al.l (1201 il l and ICaroneet^jh 
(120141 ) and make use of the perturbation method (iHolton 


I 1992 I I. We will focus on perturbations in horizontal velocity 
V and geopotential height z = ^/g as proxies for planetary 
waves. In perturbation theory, variables are divided into a 
basic state, which is assumed to be constant in time and lon¬ 
gitude, and in an eddy part, which contains the deviations 
from the average. 

We will focus in this study on temporally averaged flow 
and circulation to facilitate the interpretation of long-term 
climate states. We can therefore assume 


U = II -I- u , 
and 


( 18 ) 


z = z + z', (19) 

where the first term in each equation denotes the basic state, 
that is, the zonal mean of a given quantity, and the second 
term the eddy component, respectively. 

It will be shown that the perturbation method makes 
time averaged Rossby wave features directly visible and fa¬ 
cilitates the interpretation of the physical processes that 
produce the observed large scale circulation patterns. The 
perturbation method in combination with the meridional 
Rossby wave numbers proves particularly powerful. 

Since we focus in this study on contributions of standing 
tropical and extra tropical Rossby waves, we provide in the 
following a rough idealized sketch on expected geopoten¬ 
tial height anomalies and resulting idealized flow patterns 
(Figure [2|, both for the respective meridional Rossby wave 
numbers close to unity. 

Eddy velocities if (not shown) exhibit cyclonic move¬ 
ment around negative and anti-cyclonic movement around 
positive geopotential height anomalies (Figure O left pan¬ 
els) . Eddy wind patterns circulating around one eddy geopo¬ 
tential height minimum and maximum per hemisphere are 
so called Rossby wave gyres for zonal wavenumber n = 1. 
They are not centred at th e equator - eve n when induced by 
the tropical Rossby wave (iMatsunol 1 1 966l l . The extra tropi¬ 
cal Rossby wave gyres, however, do tend to be shifted more 
polewards compare d to their tropical counterpart (see e.g. 
lEdson et al.l (I2OIII D. In cases with meridional Rossby wave 
numbers close to unity, which implies that the planetary 
wave engulfs the whole planet, it is, however, difficult to 
distinguish between both type of waves based on their latitu¬ 
dinal location. We found in this study that the longitudinal 
location of the geopotential height anomalies are a better 
indicator. 

For tropical Rossby waves, we find Rossby wave gyres 
around negative geopotential height anomalies west from 
the substellar point and positive anomalies east from the 
substellar point (Figure [21 bottom-left). Cyclonic vorticefl 
in the observed flow patterns are thus located at longitude 
— 135° over the negative geopotential height anomaly (Fig¬ 
ure (21 bottom-right). The longitudinal day and night side 
heat gradient is predicted to trigger a tropical Kelvin wave, 
which appears in the eddy velocity field as a purely zonal 
wind strip that is ce n tred a t the equator (not shown here). 
IShowman &: Polvani! (I2OIII ') predicted that the shear be¬ 
tween the tropical Kelvin and Rossby waves produces equa¬ 
torial superrotation (Figure (21 bottom-right). 

In contrast to that, shear between the tropical Kelvin 
wave and extra tropical Rossby waves produces two west¬ 
erly jets at mid-latitude that are weaker compared to the 
purely equatorial wave interaction scenario (Figure (21 top- 
right). At the same time, the extra tropical Rossby wave 
gyres are shifted by almost 180 degrees with respect to the 
tropical Rossby wave gyres (Figure (21 top-left). As a di¬ 
rect consequence, the vortices are now located at longitude 
4-90° (Figure (21 top-right) over the now eastward shifted 
eddy geopotential height minimum. To maximize the eddy 
geopotential height signal, we will apply the perturbation 


’Cyclonic vortex’ is a type of flow that has large positive relative 
vorticity f that is defined as the vertical component of the curl 
of the horizontal velocity (^ = k X X llHoltonlll992h . 
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Figure 2. Schematic idealized view of geopotential height perturbations and flow patterns at the top of the atmosphere (p = 225 mbar) 
for the extra tropical (top) and extra tropical Rossby wave (bottom) and for meridional Rossby wave numbers close to unity, LjijR Ki 1 
and Xr/R ~ 1, respectively. Red ellipses with black continuous circumferences denote positive gepotential height anomalies, light 
blue ellipses with dashed circumferences denote negative geopotential height anomalies. Grey arrows mark the strength, direction and 
approximate latitude position of zonal wind jets, the black arrows sketch the standing planetary Rossby waves with zonal wave number 
n = 1 for each scenario. The thin black circles with arrow head mark the position of the cyclonic vortices associated with the negative 
geopotential height anomaly. 


method t o the upper atmosphere layer at p = 225 mbar, 
following tedson et al.ll2oill : 1carone et al.ll2014l) . We have 
verihed that the sponger layer at p = 100 mbar does not 
affect our results. We will also analyze the horizontal flow 
pattern at this altitude, which allows to directly tie changes 
in the Rossby waves to changes in the horizontal flow. 

We hnd in this study either climate states dominated 
by tropical Rossby waves or a mixture between tropical and 
extra tropical Rossby waves. We will show in an upcoming 
paper, however, that it is indeed possible to achieve a climate 
state strongly dominated by extra tropical Rossby waves 
with our model, if surface friction is increased. 


5 RESULTS AND DISCUSSION 

We will discuss the evolution of the horizontal flow from 
slow Prot = 100 days to very fast rotations Prot = 1 day 
and focus on the possible transition regions associated with 
phase changes in the Rossby waves excited in the planet’s 
atmosphere. The latter are associated with changes in the 
meridional Rossby wave numbers, Lr/Rp and Xr/Rr. The 
transition in meridional Rossby wave numbers are inspected 
using the perturbation method. The combination of both 
methods allows to precisely identify which type of planetary 
wave is responsible for a given climate state. 


We will show that we can identify more climate state 
transitions than previous studies due to our very thorough 
investigation of rotation period and planet size phase space. 
In particular, we identify two regions in the fast rotation 
regime that allow for the existence of multiple climate states. 
We will show that different climates states result in different 
zonal wind jet structures, longitudinal cyclonic vortex loca¬ 
tions, and different the upper atmosphere hot spot locations 
with respect to the substellar point. 


We will furthermore discuss circulation cell strength 
and number and their interaction with horizontal flow pat¬ 
terns in different climate states. We will show that fast equa¬ 
torial superrotation suppresses direct circulation and that 
large planets can exhibit fragmented circulation cells. We 
will also discuss surface temperature evolution in compari¬ 
son with the previously obtained results, in particular, with 
the diminishing strength of the tropical direct circulation 
due to strong equatorial superrotation. Finally, we will in¬ 
vestigate the vertical temperature and explain two different 
types of upper atmosphere temperature inversions based on 
the insights gathered from previous results. 
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Figure 3. Maximum zonal wind speeds in m/s on a tidally locked 
terrestrial planet for different planet radii {Rp = 1 — ‘2REarth) 
and rotation periods {Prot = 1 — 100 days), assuming identical 
thermal forcing, that is, moderate Earth-like irradiation. 


5.1 Zonal wind 

Figure[3]shows the ruEiximum zonal wind speeds with respect 
to planet rotation period and planet size. A few things are 
immediately apparent: 

• Positive values indicate that the zonal wind maxima 
are associated with westerly winds, that is, the atmosphere 
flows in general in the direction of the planet’s rotation. 

• Large planets tend to have larger wind speeds. 

• Wind speed increases for faster rotation periods start¬ 
ing with a few tens of m/s for the slowest rotation (Prot = 
100 days) and reaching velocities larger than 100 m/s for 
very fast rotation {Prot = 1 day). 

• For fast rotating planets with 1.5 Si Prot ^ 5 days, the 
wind speed maxima show a wide spread of possible values. 

• For very fast rotation, Prot = 1 day, the wind speeds ap¬ 
pear to converge for all planet sizes to wind speeds between 
275-300 m/s. 

The westerly flow ca n be understood according to 
IShowman fc Polvanil (l201ll ) by the shear between equatorial 
(eastward travelling) Kelvin waves and tropical (westward) 
travelling Rossby waves that lead to net angular momen¬ 
tum transport from high to low latitudes and thus to west¬ 
erly flow. The meridional tilt of the Rossby wave gyres, a 
clear indication of the Kelvin-Rossby wave coupling - will 
be shown later for fast rotating planets. The other zonal 
wind tendencies warrant a more detailed analysis. 


5.1.1 Dependency of zonal wind speeds on planet size 


The tendency of larger planets to have larger wind speeds 
can be understood by performing a scale analysis on the 
thermodynamic energy equation. The first law of thermo¬ 
dynamics can be written for synoptic s cales and horizontal 
temperature fluctuations jHoltorJI 19921) as 


T ( do de de\ J 
To ) ~T’ 


( 20 ) 


Table 4. Zonal wind speeds in the troposphere and expected 
wind speeds Uexp from scale analysis for Prot = 100 days and 
different planet sizes 


Rp 

\.^Earth\ 

^max 

[m/s] 

'^exp 

[m/s] 

1 

5 

5 

1.25 

11 

8 

1.45 

12 

10.5 

1.75 

18.4 

15.3 

2 

19.8 

20 


if we neglect vertical transport of energy (see also Equa¬ 
tion 0). 0o{z) is the basic state of the potential tempera¬ 
ture, 0{x,y,t) the deviation from it, where 2 is the vertical 
component in the Cartesian coordinate system, and J is the 
heating rate. Because we assume the same temperature forc¬ 
ing and atmosphere composition throughout simulations in 
this study, we can assume that dd jdt is approximately con¬ 
stant for changes with planet size. Therefore, we can derive 
the planet size dependence of the remaining wind advection 
term on the left ha nd side of equation dsni) via scale analysis 
with (lHoltonlll99j) . 

L cp 


where U is the horizontal velocity scale, and L is the horizon¬ 
tal lengt h scale. As already disc ussed in Section r2.4[ L — Rp 
(see also lShowman et al.l (I2OIII )). 9{x, y) is in the context of 
tidally locked terrestrial planets the difference between day 
and night side temperature, which we prescribed as equal for 
all planet sizes. We use Newtonian cooling for the heating 
rate, so that 



'^rad 


( 22 ) 


where A0 is the difference between the potential temper¬ 
ature 9 in the atmosphere and the prescribed equilibrium 
potential temperature 9eq. Trad is again the radiative time 
scale. If we assume that A9 hardly changes with planet sizes, 
then the heating rate J dependency on planet size is deter¬ 
mined via Trad- The radiative time scale depends inversely 
on surface gravity, where g oc Rp from Section 12.41 From 
equations (HU and ([22ll. it follows that J oc Rp. Because 
also L — Rp, it follows ultimately from equation m that 
wind speeds depend on planet size via relation: 

U oc Rp. (23) 


Therefore, we expect stronger wind speeds by a factor of 4 
for the 2 Earth radii planet compared to the 1 Earth radii 
planet. This relation, however, only holds for the free at¬ 
mosphere and if the Coriolis force is negligible, that is, for 
the atmosphere above th e planet bounda r y laye r and slow 
rotating planets (see, e.g. IShowman et al.l (1201 ll )). 

For Prot = 100 day, the maximum zonal wind speed 
in the troposphere of the IRparth planet is Umax = 5 m/s 
and for the 2REarth planet it is Umax = 19.8 m/s. This 
is very close to the 1:4 relation expected from the scale 
analysis. The R% increase in zonal wind speeds is, indeed, 
maintained for the investigated planet sizes within a few 
m/s if Umax ~ 5 m/s is taken as a basis value from the 
Rp = IRparth simulation (Table |4]). This is a good agree- 
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Latitude [degrees] 


Figure 4. Zonal wind speeds in m/s in our nominal climate model for Prot = Porb = 1 ds-y- Minimum and maximum wind speeds are 
-20 to 280 m/s, respectively. Contour interval is 10 m/s. 


ment between the scale analysis and the simulation results, 
if we assume that Umax oc U. For faster rotation, however, 
pure advection from the hot day to the cold night side is 
no longer a valid assumption as the Coriolis force becomes 
a more dominant factor in shaping the flow. While bigger 
planets still tend to have higher wind speeds for faster ro¬ 
tation, the differences are much smaller than expected from 
U ~ B?p, and even appear to vanish for very fast rotation. 


5.1.2 Dependency of zonal wind speeds on planet rotation 


The increase of zonal wind speed with faster rotation is more 
complex and warrants an in-depth investigation. Generally, 
we see an increase of zonal wind speeds with faster rotation 
that reach maximum values of 300 m/s for Prot = 1 day. 
The high wind speeds of the very fast rotating planets can 
be directly attributed to equatorial super-rotation. Figure |4] 
shows exemplary zonal mean of zonal winds for Prot = 1 day 
and different planet sizes. 

Equatorial superrotation is a dynamical feature that 
frequently and consistently shows up in the atmosphere 
modelling of tidally l ocked hot Jupiters and Neptunes (e.g. 


modelling or tidally l ocked not Jupiters and INeptunes fe.g. 
iKataria et ^ (|201^), Lewis et al.| J 2 OI 0 I) ) and for hot Super- 
Earths fe.g. lKataria et aTT i 2014ll . lZalucha et al.l (l2013l ll. In¬ 
terestingly, atmosphere models of fast rotating tidally locked 
terrestrial exoplanets with habitable surface temperatures 
and the same nitrogen-dominated atmosphere composition 
are ambiguous with regard t o equa torial superrotation. The 
simulations of lEdson et al.l ll201lll also show zonal wind 


speed increase with faster rotation, but hardly exceed ve¬ 
locities of 100 m/s. Their simulations show two jets at high 
latitudes v = ±60° for fast rotation {Prot ^ 4 days) in¬ 
stead of one equatorial jet. Even for Prot = 1 day, their 
simulations still have high latitude jets in addition to an 
equatorial superrotating jet w i th tw o maxima located at 
u = ±10°. iMerlis fc Schneider! (l2010ll performed a simula¬ 
tion for a tidally locked ocean planet with Prot = 1 day 
that exhibited two zonal jets at high latitude {v = 60°) 
and no equatorial superrotating jet. Both studies assumed 
an Earth-size planet. In contrast to that, all our simulations 
show efficient equatorial superrotation across all planet sizes. 

IShowman fc Polvanil (1201 il l attributed the emergence 
of equatorial superrotation and the associated fast wind 
speeds to the coupling between standing equatorial Kelvin 
and Rossby waves, provided the tropical Rossby radius of 
deformation Xr (Equation I17II is equal to or smaller than 
the planetary radius. 

lEdson et al.l (1201 ill . on the other hand, speculated that 
the extra tropical Rossby radius of deformation Lr (Equa¬ 
tion 0 may be an important factor in shaping atmosphere 
dynamics in their simulations. They identified an abrupt 
transition in large scale circulation and zonal wind speeds 
at Prot = 4.2 days, which coincides with Lr/Rr ~ 1 . Thus , 
the phase transition occurs according to lEdson et al.l (1201 ill 
because the extra tropical Rossby waves starts to ’fit’ on the 
planet. It can be seen from their Table 2 that this phase tran¬ 
sition coincides with another transition: Xr/Rr « 0.5. The 
latter may hint to an involvement of the equatorial Rossby 
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Figure 5. Left panel: Equatorial Rossby radius of deformation Xr over planet radius Rp. Right panel: Mid-latitude Rossby radius 
of deformation evaluated at 30° latitude over planet radius Rp. Both properties are shown for different planet sizes {Rp = 1 — “^REarth) 
and rotation periods {Prot = 1 — 100 days). 


wave as well. For rotations slower than the rot ation period 
at the transition, that is, for Prot = 5 days, lEdson et al.l 
ll201l]l explicitly shows the presence of equatorial Rossby 
gyres in the eddy geopotential heights (their Figure 5g) and 
broad equatorial zonal wind jet feature (their Figure 4g). Af¬ 
ter transition {Prot = 4 days), the zonal wind jet and gyre 
structure abruptly changes to a high latitude configuration 
centred at iz = 60° (their Figures 4e and 5e). 

In the light of the differences in wind jet structure 
for fast rotating planets from different climate models, we 
will study climate dynamics for different rotation periods 
in what follows. We will focus on standing Rossby waves 
and will show that we can link climate state transitions 
to either tropical or extra tropical Rossby wave transitions. 
The transitions can be inferred from the meridional Rossby 
wave numbers and further investigated with the perturba¬ 
tion method by means of the eddy wind and eddy geopo¬ 
tential structure at p = 225 mbar. We will, furthermore, in¬ 
vestigate how climate state transitions affect horizontal flow 
structure, where we also focus on the upper atmosphere at 
p — 225 mbar. 


5.2 Phase transition in horizontal flow with 
meridional Rossby wave number 

As already pointed out bv lEdson et al.l (I 2 OIII I. it is not pos¬ 
sible to exactly predict the behaviour of Ljt and An for a 
3D planet atmosphere as both Rossby radii of deformation 
depend also on the vertical temperature structure via dO/dp 
and on the surface temperatures via the scale height H. Both 
dynamical properties change as the atmosphere dynamics 
changes, as will be shown later. The meridional Rossby wave 
numbers Lr/Rp and Xp/Rp not only decrease with faster 
rotation, they also depend on planet size (Figure [5] and Ta¬ 
ble O. 

Differences in surface te mperature and vertic al struc¬ 
ture, as already pointed out in ICarone et al.l ll2014h . explain 
why our Lp/Rp « 1 transition appears at slower rota¬ 
tion, Prot = 6 — 7 days, f or Rp = ^REarth compared to 
the model of lEdson et al.l ll201ljf that shows transition at 
Prot = 4.2 days (Figure [5] and Table [5|. 

In the following, we will examine horizontal flow pat- 


Table 5. Rossby wave number transition for different planet sizes 
in our nominal model. 


Rp 

Prot for 

Rp ~ 

0.5 

Prot for 
Lr 1 

Rp~ 

Prot for 

Rp~ 

Prot for 

>'R 

Rp ^ 

0.5 

\_REarth\ 

[days] 

[days] 

[days] 

[days] 

1 

3-4 

6-7 

25 

5-6 

1.25 

4 

8-9 

30 

7 

1.45 

5 

10 

34 

8 

1.75 

5-6 

11-12 

35-40 

9 

2 

5-6 

12 

40-45 

10-11 


terns at four transition regions in meridional Rossby wave 
numbers (Tabled, where we will show that the second and 
third transition almost coincide and can thus not be inves¬ 
tigated separately: 

• \p/Rp «1 

• Lp/Rp «1 

• Xp/Rp ~ 0.5 

• Lp/Rp ~ 0.5 

5.2.1 Rossby wave transition region: Xp/Rp ~ 1 

The transition region where the meridional extent of the 
tropical Rossby wave starts to ’fit’ on the planet can be 
found between 25 days for an Earth-size and 45 days for a 
planet twice the size of Earth (Table [5|). When comparing 
with Figure (3] this region coincides also with an increase 
in zonal wind speed acceleration with shorter rotation peri¬ 
ods: The Xp/Rp > 1 region with Prot = 45 — 100 days, is 
a region of low zonal wind speeds between 10-30 m/s. The 
Xp/Rp ^ 1 region with Prot = 20 — 35 days, shows an in¬ 
crease in zonal wind speeds to velocities of up to 50-75 m/s, 
depending on the planet size. This acceleration is in line with 
the interpretation that the equatorial Rossby waves form 
standing waves and transport angular momentum more effi - 
ciently from high to low latitudes. Also lLeconte et ^ (l2013ll 
pointed out that their 3D simulations show direct circulation 
states for Xp/Rp >> 1 . 

Inspection of the eddy geopotential height z' and eddy 
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Figure 6. Left panel: Eddy geopotential height in [m] and eddy horizontal winds in [m/s] at p = 225 mbar. Right panel: Zonal mean 
of zonal winds in [m/sj. Both are shown for an Earth-size planet (Rp = IREarth) rotation periods at the Rossby wave transition 
region Xp/Rp ~ 1 {Prot = 32,35, and 45 days for top to bottom rows). The longest wind vectors are from top to bottom: 27.6, 24.3, 
and 32.2 m/s, respectively. Contour intervals are 10 m for the eddy geopotential and 1 m/s for zonal wind speeds, respectively. 



Figure 7. Temperatures (contour interval 10 K) and streamlines 
of the horizontal flow, averaged over 1000 d, for pressure level p = 
225 mbar for an Earth-size planet {Rp = ^Rparth) rotation 
periods at the Rossby wave transition region Xp/Rp ~ 1 (from 
top to bottom: Prot = 32,35, and 45 days). 


horizontal winds (see Section 0 at the top of the atmo¬ 
sphere p = 225 mbar shows indeed in this planet rotation 
region a gradual change from a climate state dominated by 
direct circulation to one where a standing tropical Rossby 
wave appears and starts to contribute to angular momentum 
transport - hence the zonal wind speed increase (Figure 0 
left panel). Direct circulation is identihable by divergent flow 
and a maximum in z' at the substellar point, which marks 
the top of the direct circulation cell with horizontal trans¬ 
port from the hot substellar point towards the cool night 
side aloft and upwelling flow at the hot substellar point be¬ 
low. The standing tropical Rossby wave is identihable by the 
Rossby wave gyres with negative eddy geopotential height 
east of the substellar point and positive west of it (See Fig¬ 
ure 0). 

_ One of the climate states discussed in ICarone et al.l 

(|2014| . the Super Earth planet with Rp — 1.4:5REarth, 
g = 14.3 m/s^ and Prot = 36.5 days resides in the transition 
regime. It was likewise shown to exhibit a mixture between 
superrotating flow and divergent flow at the top of the tro¬ 
posphere. It was then tentatively assumed that there is a 
connection between cyclonic vortices, the standing tropical 
Rossby wave and the zonal wind structure. This connection 
can be conhrmed in this study: Although, weak cyclonic vor¬ 
tices are present also for slow rotations (e.g. Prot=45 days), 
the cyclonic vortices gain in strength and become locations 
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associated with upper atmosphere temperature maxima at 
the same time as the planetary wave starts to form a stand¬ 
ing wa ve (Figure El right p anel and Figure [7]). 

In ICarone et alJ (l2014l h we speculated that the temper¬ 
ature maxima associated with the vortices (Figure[7]) arouse 
from adiabatic heating of air masses streaming from the 
substellar point downwards into the vortices; The substellar 
point is a region of particularly high geopotential height. At 
the same time, the air emerging at the upper troposphere 
over this point is particularly cool, which can be explained 
by adiabatic cooling during the upwelling process. These 
cool air masses ’fall’ downwards towards the vortices, which 
are cold features i n the mid-atmosphere at p = 525 mbar 
(see lCarone et ahl ll20l4l ') and show up as negative anoma¬ 
lies in the geopotential height held. Consequently, the air is 
adiabatically heated over the vortices. We will show in Sec- 
tion l7.3l that the temperature maxima over the vortices cause 
temperature inversions in the vertical temperatu re struc¬ 
ture, as already pointed out in ICarone et al.[ (l2014l l. 

The transition to a standing tropical Rossby wave 
regime is also accompanied with a transition in wind struc¬ 
ture: From a weak equatorial superrotation towards a broad 
zonal wind structure with two zonal wind maxima at mid¬ 
latitude and mid-troposphere (p = 525 mbar) (Figure 
right panel). Apparently, angular momentum is indeed more 
efHciently transported from the poles towards lower latitudes 
with transition to Xr/Rp ^ 1, as expected from shear be¬ 
tween tropical Rossby waves and Kelvin waves. The latter 
are visible in the eddy velocity wind held as very conhned 
equatorial stripe (between p = ±10°) of purely zonal eddy 
winds, except over the substellar point where divergent up- 
welling how superposes the Kelvin wave structure (Figure|6l 
left panel). 

According to Figure [O] the emergence of the stand¬ 
ing tropical Rossby wave can be found by visual inspec¬ 
tion at Prot = 35 days instead of Prot = 25 days based on 
Xr/Rp (Table (Sjl. Prot = 35 days corresponds, however, to 
a Xr/Rp « 1.15, which is only slightly larger than unity. 
Thus, we conhrm that Rossby wave numbers give a good es¬ 
timate but albeit not an exact guideline to identify changes 
in climate states. 

Comparison between the eddy geopotential height for 
Rp = IRparth and Prot = 35 days (Figure [6]), with the 
eddy geopotential height for Rp = 1.4:5REarth and P rot = 
36.5 days (Figure 7, right panel in ICarone et al.l (|2014l ')') con- 
hrms the planet size dependency. The Rossby wave gyres 
are more apparent for the Super-Earth-size planet than 
for the Earth-size planet, even thoug h the rotation of the 
tidally locked Super-Earth discussed in I Carone et al.1 (l2014l l 
is slightly slower. This is in line with the prediction that the 
rotation period for which standing equatorial Rossby waves 
emerge is shifted towards slower rotations with increasing 
planet size. 

5.2.2 Rossby wave transition region: Lr/Rp « 1 and 
Xr/Rp « 0.5 

This is the regime that lEdson et all (l201ll l associated with 
multiple solutions for the climate state and with abrupt 
transi tion in the zonal wind speeds. Like in lEdson et al.l 
(|2nilll . the tropical and extra tropical Rossby wave transi¬ 
tion appear to almost coincide for all planet sizes (Table [5]). 


Furthermore, another possible phase transition for the ex¬ 
tra tropical Rossby wave with Lr/Rp ~ 0.5 is not very 
far away for an Earth size planet at Prot = 3 — 4 days. 
For the largest planet in our sample {Rp = 2REarth), the 
Lr/Rp « 0.5 transition lies at Prot =5-6 days and thus far¬ 
ther away from the Lr/Rp « 1 transition at Prot = 12 days 
and the Xr/Rp « 0.5 transition at Prot = 10—11 days, re¬ 
spectively. We will thus focus in the following discussion on 
atmosphere dynamics on the largest planet in our sample. 

Inspection of zonal wind speed maxima (Figure [3)1 
shows that while the velocities still generally increase with 
faster rotation, the increase becomes less steep between 
9 ^ Prot ^ 20 days. As soon as the equatorial Rossby wave 
transits to Xr/Rp < 0.5 , the wind speeds rise again steeply. 
We do not End, however, an abrupt transition in wind speeds 
as reported bv lEdson et al.1 (|201ll L This diflerence in zonal 
wind evolution will become clearer in the next section, as 
we explore why our mo del goes into full equ atorial superro¬ 
tation and the model of lEdson et al.l (1201 il l does not. 

Figures [8] and |9] show indeed a change in climate state 
near the Xr/Rp « 0.5 transition that affects the upper 
branch of the direct circulation cell. Upwelling from the 
substellar point can again be recognized as a cold spot at 
p = 225 mba]0 (Figure [O] bottom panel) and as a positive 
geopotential height anomaly with outstreaming flow in the 
eddy geopotential height and winds (Figure [8l left-bottom 
panel) for Prot = 16 days. After Xr/Rp « 0.5, upwelling and 
divergence weaken and disappear as the Rossby wave gyres 
become even stronger (Figure lU top panel and Figure (Sj 
left-top panel). At the same time, the zonal wind jet max¬ 
ima converge towards the equator and increase in strength 
(FigurelH right panel). It will be shown in Section |6] that di¬ 
rect circulation and strong equatorial superrotation, indeed, 
don’t mix very well: If direct circulation is strong, superro¬ 
tation is diminished and vice versa. 

With the disappearance of the divergence, the tem¬ 
perature of the cyclonic vortices at p = 225 mbar - that 
have been pushed towards the poles with faster rotation - 
changes drastically. For Prot = 16 days, the vortices are still 
relatively warm, for Prot = 9 days, however, they are in¬ 
stead associated with very cool temperatures even at this 
low pressure levels. This result can again be explained in 
the context of divergent flow emerging from the substellar 
point ’falling down the slopes’ of the warm - and thus verti¬ 
cally extended - equatorial superrotating region. Upwelling 
flow that falls into the vortices with low geopotential height 
at the top of the troposphere becomes, thus, adiabatically 
heated (see previous Section): provided the divergent flow 
can reach the increasingly pole-wards shifted vortices and 
provided there is divergent flow to heat up in the first place. 
The removal of divergent flow with faster rotation also ex¬ 
plains why the horizontal temperature gradient increases be¬ 
tween Prot = 16 days and Prot = 10 days at the top of the 
atmosphere (Figure [9] and also Figure [3 for comparison). 
We will take up the discussion of upper atmosphere thermal 
inversion again in Section 17.31 

The atmosphere dynamics of our Prot = 9 days sim¬ 
ulation for Rp = 2REarth is, furthermore, very similar to 

® As noted in I Carone et al.l ||2014^ , the cool temperatures arise 
due to adiabatic cooling during uplifting. 
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Figure 8. Eddy geopotential height in [m] and eddy horizontal wind in [m/s] at p = 225 mbar (left panel) and zonal mean of zonal 
winds in [m/s] (right panel) for a large Super-Earth planet (Rp = ^Rparth) ^md rotation periods at the Rossby wave transition region 
LpilRp fa 1 and XpjRp fa 0.5 {Prot = 9,13, and 16 days). The longest wind vectors are from top to bottom: 62.6, 59.0, and 61.1 m/s, 
respectively. Contour intervals are 50 m for the eddy geopotential and 2 m/s for zonal wind speeds, respectively. 



Figure 9. Temperatures and streamlines of the horizontal flow, 
averaged over 1000 d, for pressure level p = 225 mbar for a large 
Super-Earth planet {Rp = “iREarth) ^md rotation periods at the 
Rossby wave transition region Lp/Rp fa 1 and Xp/Rp fa 0.5 
(from top to bottom: Prot = 9,13, and 16 days). Contour intervals 
are 2 K for Prot = 16 days and 5 K for Prot = 9, and 13 days. 


the Prot = 10 days simulation in ICarone et al.l (|2014l ) for 
Rp = lAbREarth and the Prot = 5 days simulation in 
lEdson et al.l 1 201ll ). All these cases lie near their respec¬ 
tive \p/Rp = 0.5 transition regime, which shows again how 
valuable the meridional Rossby wave numbers are as a di¬ 
agnostic tool as they allow to compare results from different 
climate models with each other. 

The increase in zonal wind speeds, the evolution of the 
zonal wind jet towards the equator and the deepening of 
the tropical Rossby wave gyres in the eddy wind and eddy 
geopotential height field show that the standing tropical 
Rossby waves are still the driving mechanism in this climate 
state regime. This is remarkable because the extra tropical 
meridional Rossby wave number {Lp/Rp ^ 1) shows that 
also the formation of standing extra tropical Rossby waves 
is, in principle, possible. The latter apparently don’t yet play 
a role - at least in our nominal model. We will show in an 
upcoming paper, that an increase in surface friction trig¬ 
gers indeed the formation of standing extra tropical Rossby 
waves for Lp/Rp $( 1. 


5.2.3 Rossby wave transition region: Lp/Rp « 0.5 

The next climate state transition occurs for a large terres¬ 
trial planet with Rp — 2REarth at Prot = 5 — 6 days. Indeed, 
the eddy geopotential height and eddy winds as well as the 
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Longitude [degrees] Longitude [degrees] 

Figure 10. Eddy geopotential height in [m] and eddy horizontal wind in [m/s] at p = 225 mbar (right panel) and zonal mean of zonal 
winds in [m/s] (left panel) for a large Super-Earth planet {Rp = 2Rparth) rotation periods at the Rossby wave transition region 
LrIRp ~ 0.5 {Prot = 3, 4, and 6 days). The longest wind vectors are from top to bottom: 61.5, 54.7, and 63.5 m/s, respectively. Contour 
intervals are 25 m for the eddy geopotential and 5 m/s for zonal wind speeds, respectively. 
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Figure 11. Temperatures and streamlines of the horizontal flow, 
averaged over 1000 d, for pressure level p = 225 mbar, Rp = 
2Rp^^i^ and rotation periods at the Rossby wave transition re¬ 
gion Lp/Rp fa 0.5 (from top to bottom: Prot = 3, 4, and 6 days). 
Contour intervals for temperatures are 5 K. 


zonal wind structure show again a change in dynamics be¬ 
tween Prot = 6 days and Prot = 4 days (Figures [10] and 

HD). 

Extratropical Rossby wave gyres appear besides the 
tropical Rossby wave gyres for Prot = 4 days. These can be 
identified by the apperance of additional geopotential height 
anomalies at high latitudes that are furthermore shifted in 
longitude by almost 180° with respect to the tropical Rossby 
wave gyres (See also Figure [2|). That is, the extra tropical 
positive z' structure is located west from the substellar point 
and the negative east from the substellar point, respectively 
(Figure (To] middle panel). 

The extratropical standing Rossby waves lead to the 
formation of additional high latitude wind jets at iz = 60° 
(Figure flUl and im middle panel). Because the equatorial and 
the extratropical Rossby wave have both meridional extent 
of less than half the planetary radius, both waves can now 
form standing planetary waves. The existence of two plan¬ 
etary waves can also be inferred from the horizontal flow 
streamlines at the upper troposphere (Figure 1111 top and 
middle panel). 

For even faster rotation, Prot = 3 day (Figure ITO] top 
panel)), the equatorial standing Rossby waves become again 
dominant, although remnants of the extratropical Rossby 
wave remain visible in the eddy geopotential wind field: We 
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Figure 12. Eddy geopotential height in [m] and eddy horizontal winds in [m/s] at p = 225 mbar (top panels) and zonal mean of zonal 
winds in [m/s] (bottom panels) for Rp = IRparth (left) and 1.25REarth (right) and rotation period Prot = 3 day. The longest wind 
vectors are from left to right: 39.6, and 37.4 m/s, respectively. Contour intervals are 50 m for the eddy geopotential and 1 m/s for zonal 
wind speeds, respectively. 



Longitude [degrees] Longitude [degrees] 

Figure 13. Temperatures and streamlines of the horizontal flow, averaged over 1000 d, for pressure level p = 225 mbar for Rp = IRparth 
(left) and l.2bREarth (right) and rotation period Prot = 3 days. Contour intervals for temperatures are 2 K. 


still find positive high latitude geopotential anomalies be¬ 
tween longitudes —180° to -135° and negative geopotential 
anomalies east from the substellar point between longitudes 
-1-60° to -1-135° that appear to supersede the negative geopo¬ 
tential anomaly associated with the tropical Rossby wave 
(between longitudes —45 to +135°). We also find in the hor¬ 
izontal flow that the longitudinal location of the vortices is 
flipped with respect to the substellar point (Compare Fig¬ 
ure EH top and middle panel with Figure O): They are still 
located at —135° for Prot = 6 days and are shifted by the 
appearance of the extra tropical Rossby eave to longitudes 
between +90° to +135° for Prot = 4 and 3 days, respectively. 

Strikingly, the form and meridional tilt of the Rossby 
wave gyres bet ween —45° and 135° longit ude are exactly 
as predicted bv IShowman fc Polvanil (l201ll) (Compare their 
schematic Figure 1 with Figure EOl top-left panel). As pre¬ 
dicted, these tilted standing Rossby waves lead to a very ef¬ 
ficient spin-up of equatorial superrotation for Prot ^ 3 days 
(E.g., Figure[3]). It will also be shown in the subsequent sec¬ 
tion that remnants of a weak standing extra tropical Rossby 
wave can be found, even for very fast rotating planets with 
strong equatorial superrotation. 

Interestingly, the simulations of lEdson et al.l (l201ll ) 


don’t go into a dynamics regime that is fully dominated by 
standing tropical Rossby waves for Prot ^ 4 days (See their 
Figures 4a,c,e and Figures 5a,c,e). Instead, their Prot = 4 
and 3 days simulation are dominated apparently by standing 
extra tropical Rossby waves and correspondingly by 60° high 
latitude jets and a pair of very weak jets near the equator. 
Their Prot = 1 day simulation shows a mixture of extra trop¬ 
ical and tropical Rossby wave gyres. Again, these Rossby 
waves drive two zonal jets per hemi sphere: one a t the e quator 
and one at high latitudes. Because lEdson et al.l (1201 ill never 
find a purely superrotating state, the a uthors report maxi¬ 
mum wind speeds of at most 111 m/s. iMerlis fc Schneider! 
ll2010l) have no equatorial superrotation in their Prot = 1 day 
simulation and find consequently even slower wind speeds 
of about 50 m/s. In contrast to that, our simulations reach 
300 m/s for Prot = 1 day and all planet sizes (Figure [4]). 

As noted in the previous section, lEdson et al.l (l201l[l 
identified an abrupt transition in zonal wind speeds and at 
least two dynamical state solutions between Prot =3-4 days 
in their dry simulation. The authors tentatively linked the 
rotation region of multiple solutions to a Lp/Rp = 1 transi¬ 
tion. Now we can explain why we didn’t see such an abrupt 
transition in zonal wind speeds and also go into a different 
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climate state for faster rotations: The extra tropical Rossby 
w ave is stronger than tropical Rossby waves in the model 
of lEdson et alJ (|201ll '). In our model, the equatorial Rossby 
wave is dominant. Th e dominance of ext ra tropical Rossby 
waves in the study of lEdson et al.l fSoTJ) also explains the 
abrupt transition in zonal wind speeds and structure - from 
equatorial to high latitude - as soon as the meridional ex¬ 
tent of the extratropical Rossby wave ’fits’ on the planet 
and allows to form standing waves for Prot ~ 4 days. In our 
simulation, equatorial Rossby waves remain dominant after 
transition. The extra tropical Rossby waves are suppressed 
until both X/Rp and Lp/Rp become both smaller than 0.5 
and both waves can ’fit’ on the planet. However, even in this 
case the tropical Rossby wave is dominant most of the time. 

Thus, we conclude that there are two possible climate 
states in the 0.5 ^ Lp/Rp ^ 1 regime: One that is domi¬ 
nated by standing tropical planetary waves with equatorial 
fast zonal jets - as in our simulation - and one that is domi¬ 
nated by standing extra tropical planetary waves wi th slower 
high latitude jets - as shown in lEdson et al.l 1 I 2 OIIII. We as¬ 
sume that the multiple equilibria mentioned bv lEdson et al.l 
(120111) refer to these two states. 

The picture becomes even more complicated for 
Lp/Rp ^ 0.5, where we conclude that there are now at 
least three dynamical states possible: a state dominated by 
tropical Rossby waves, a state dominated by extra tropical 
Rossby waves, and a mixture of both type of waves. As soon 
as there is an extra tropical wave component, the zonal wind 
is slower compared to the purely tropical Rossby wave state. 

The existence of different climate states makes the com¬ 
parison of zonal wind speeds in the Lp/Rp ^ 0.5 rotation 
regime difficult and also explains the previously identified 
spread in zonal wind maxima (Figure [S]): Larger planets 
don’t always have higher wind speeds than smaller plan¬ 
ets and there appear to be ’jumps’ in the evolution of zonal 
wind velocity with faster rotation. Only climates with strong 
equatorial Rossby waves experience very fast winds as the 
angular momentum is transported from high latitudes all 
the way towards the equator. Extra tropical waves divert 
at least part of the angular momentum towards 60° lati¬ 
tude. Thus, climate states with high latitude jets are always 
slower than fully formed equatorial jets. Because small plan¬ 
ets tend to have smaller wind speeds at the same rotation 
period than bigger planets. Earth-size planets that exhibit a 
mixture of tropical and extra tropical Rossby waves should 
show particular low wind speeds. Indeed, the slowest zonal 
winds in the Lp/Rp « 0.5-regime with Umax ~ 50 m/s can 
be found for the smallest planets in our study, Rp = 1 and 
Rp = 1.25REarth, where both simulations show a mixture 
between equatorial and extra tropical wave dynamics (Com¬ 
pared Figures \TI\ and [TSl with Figure [SJ. 

5 . 2.4 Very fast rotation: Prot = 1 day. 

As already shown in Figure [d] our simulations show for all 
planet sizes strong equatorial superrotation. Simulations of 
tidally locked large exoplanets, that is, hot Jupiters (e.g. 
iKataria et al.l ((201^) and even larger Super-Earth planets 
than investiga t ed in this study (e.g. for Rp — 2.75REarth, 
IZalucha et al.l (l2013l )) show anti-rotating, that is, easterly 
jets at higher latitudes. Indeed, the two largest planets in 
this study, Rp = 1.75REarth and Rp = 2REarth show signs 


Table 6. Latitudinal extent of Rossby radius of deformation for 
Prot = 1 day in comparison with real width of equatorial jet. 


Rp 

[^Earth] 

[degree] 

Lr 

[degree] 

eq. jet width 
[degree] 

1 

± 20.3 

± 16.7 

± 35-45 

1.25 

± 17.3 

± 11.2 

± 30-40 

1.45 

± 15.3 

± 9.3 

± 30 

1.75 

± 14.3 

± 8.0 

Ri ± 25 

2 

± 13.7 

±7.2 

Ri ± 20 


of anti-rotating jet formation at the surface that become 

more pronounced with planet size. _ 

According to IShowman fc Polvanil 1 I 2 OIII) . the latitudi¬ 
nal width of an equatorial superrotating jet should corre¬ 
spond to the equatorial Rossby radius of deformation. This 
is ±20 degree for an Earth-sized planet and ±13 degree for 
a planet twice the size that both rotate with Prot = 1 day 
(FigureO left panel and Table[6ll. On the one hand, there is 
indeed a stronger confinement of the equatorial jet apparent 
for larger planets (Figure |4] and Figure [Till . On the other 
hand, the jet widths are consistently broader than expected 
if the equatorial jet is formed by a standing tropical Rossby 
wave alone: The jet width at the top of the atmosphere ex¬ 
hibit a latitudinal extent of ±45 and ±20 degree for Rp — 1 
and Rp = 2RE arth, respectively (fro m visual inspection of 
Figure fT^ . Also lZalucha et aH ll2013h noted that the equato- 
rial jet width is large r in th eir simulations than predicted by 
IShowman fc Polvanil l|201ll ). Inspection of the eddy geopo¬ 
tential height shows that there are, indeed, still components 
of extra tropical wave activity - although equatorial Rossby 
wave gyres dominate the flow dynamics (Figure 1151 white 
arrows). Apparently, the positive extra tropical and equa¬ 
torial Rossby wave gyres merge around the substellar point 
for Rp = 1.25 —2REarth. The merging of equatorial and ex¬ 
tratropical Rossby waves in the strongly superrotating case 
explains why the equatorial jet width is broader than can 
be accounted for by just assuming equatorial Rossby radius 
of deformation. If we assume instead of PXp, ± {Xp ± Lp) 
as the latitudinal extent of the equatorial jet, we reach a 
better agreement (Table 1^. Thus, it appears as if the west¬ 
erly equatorial jet consists in fact of three westerly jets - 
one at the equator and one at each poleward flank, respec¬ 
tively. This interpretation is further supported by closer in¬ 
spection of the zonal wind structure, in particular for the 
largest planets investigated in this study (Figure IT}: There 
is a ’kernel’ of high wind speeds between u =100-300 m/s at 
the equator with very confined latitudinal extent. This ker¬ 
nel is flanked at each side by a broader region with slower 
winds (u = 0 — 100 m/s). 


6 CIRCULATION WITH RESPECT TO 
PLANET SIZE AND ROTATION PERIOD 

In the context of habitability, circulation is a dynamical fea¬ 
tures of uttermost importance: Circulation determines re¬ 
gions of cloud formation and of high and low precipitation. 
First insights into the change of c irculation regime wi th ro¬ 
tation period was already found in I Carone et al.l ll2014h . The 
study showed that a slow rotating tidally locked Super-Earth 
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Figure 14. Temperatures and streamlines of the horizontal flow, averaged over 1000 d, for pressure level p = 225 mbar for Rp = 
1 — 2Rparth rotation period Prot = 1 day. Contour intervals for temperatures are 5 K. The black arrows denote the hot temperature 
chevron shape. 
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Figure 15. Eddy geopotential height in [m] and eddy horizontal wind in [m/s] at p = 225 mbar for Rp = 1 — 2Rparth rotation 
period Prot = 1 day. The longest wind vectors are from left to right and from top to bottom: 110.1, 102.0, 90.0, 95.2 and 108.4 m/s, 
respectively. Contour intervals are 50 m. White arrows mark extratropical eddy geopotential height anomalies. 
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Figure 16. Left panel: Maximum strength of the clockwise equatorial circulation cell on the North hemisphere Right panel: Maximum 
strength of a counter-clockwise secondary circulation cell on the North hemisphere. The strength of the equatorial circulation cell is shown 
for the full rotation regime (Prot = 1 — 100 days), the strength of the secondary cell is shown for short rotation periods {Prot = 1 — 10 days) 
of a tidally locked terrestrial planet for different planet radii {Rp = 1 — 2Rparth)y assuming moderate Earth-like irradiation. 


planet with Prot = 36.5 days exhibits two strong direct cir¬ 
culation cells - one for each hemisphere - with an uprising 
branch at the substellar point that transports hot air to¬ 
wards the night side across the poles and terminator. This 
finding is in agreement with the majority of circulation stud ¬ 
ies on terrestrial planets (e .g. iMerlis fc Schneider! ll2010l) . 
iNavarra fc Boccalettil ll2002l) ). For Prot ~ 10 days, the cir¬ 
culation was found to be more complex: The atmosphere 
still exhibits two large global circulation cells but also con¬ 
tains two smaller embedded cou nter-circulat i ng cel ls in the 
mid-atmosphere at p =500 mbar. ljoshi et aP (Il997l ) showed 
a very similar cir culation structure for Prot ~ 16 days. 
lEdson et al.l (l201ll) . on the other hand, did not report this 
circulation state. Interestingly, even between complex Earth 
models, there is a disagreement of the circulation stat e for 
fast rotating tidally locked planets. lEdson et al.l (1201 ll) find 
three circulation cells per hemisphere for Prot = 1 day, where 
the equatorial dir ect cell is much stronger th an the other two 
circulation cells. iMerlis fc Schneider! (120101 ) report for the 
same rotation period also three circulation cells per hemi¬ 
sphere. But the polar circulation cells are the strongest in 
their simulations. 

In the following, we will precisely identify transition be¬ 
tween different circulation states with respect to rotation pe¬ 
riod and planet size. Furthermore, we will discuss possible 
implication for habitability. Henceforth, we call the circula¬ 
tion state with one unperturbed global circulation cell per 
hemisphere that is found for slow rotators state 0. This study 
allows to exactly locate the rotation periods where the em¬ 
bedded reverse circulation cells appear {state 1 ) and also to 
check if and when two or more large scale circulation cells 
emerge. We call a circulation state with two fully formed cir¬ 
culation cells per hemisphere state 2 and circulation states 
with three and more circulation cells per hemisphere state 
3. The rotation periods of regime transitions are shown in 
Table [3 The circulation states and directions of circulation 
are shown for an Earth-size planet in Figure [TT] 

The appearance of an embedded circulation cell 
(state 1) can already be identified by monitoring the maxi¬ 


Table 7. Circulation state transitions for different planet sizes in 
our nominal model. 



state 0 

state 1 

state 2 

state 3 

Rp 

Prot 

Prot 

Prot 

Prot 

[^Earth] 

[days] 

[days] 

[days] 

[days] 

1 

100 - 22 

20 - 13 

12 - 1.5 

1 

1.25 

100 - 24 

22 -12 

11 - 2 

1.5 - 1 

1.45 

100 - 26 

24 - 14 

13 - 3 

2 - 1 

1.75 

100 - 30 

28 - 10 

9 - 4 

3 - 1 

2 

100 - 34 

32 - 10 

9 - 5 

4 - 2 


mum and minimum of the meridional mass transport stream 
function Ik. This is defined as 


Vk = 


27vPp 


vdp , 


(24) 


where v is latitude and v is the zonal and temporal mean of 
the meridional velocity component u at a given latitudfl Ik 
is positive for clockwise circulation and negative for counter¬ 
clockwise circulation. 

The transition from state 0 to state 1 corresponds to 
a steep decrease in the maximum meridional mass stream 
function '^max (Figure [161 Isft panel). It sets in after the 
tropical Rossby wave transition, Xp/Rp = 1 (Tableland 
Figure [S|. Apparently, the embedded cells appear when the 
tropical Rossby wave can form a standing wave. The con¬ 
nection between tropical Rossby waves and embedded cell 
become s even strong e r whe n taking into account the find¬ 
ings in ICarone et al. (l2014l ) : The embedded reverse circu¬ 
lation cells in the meridional mass transport function ap¬ 
pear to be latitudinal cross sections of Walker-like circula¬ 
tion cells. They were shown to transport air masses primar¬ 
ily along the equator. Indicati ons of this longitud inal circu¬ 
lation could also be found in Ijoshi et al.l (1 19971 ) that also 


® Note that the horizontal wind velocity is v = (u,v), where u is 
the zonal and v the meridional component. 
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Figure 17. Meridional mass flux stream function ^ in units of 10® kg/s for Rp = IREarth showing the four possible circulation states 
and the general reduction of circulation strength with faster rotation. Note that positive values indicate clockwise and negative values 
counterclockwise circulation. Contour intervals are 5 X 10® kg/s for Prot = 1 day and 10 X 10® kg/s elsewhere. White dashed arrows 
mark the direction of circulation. 


report direct meridional circulation cells with embedded re¬ 
verse cells. A climate state with tropical Rossby waves and 
Walker-circulation can be found on Earth in the tropical 
regions and is driv en by longitudinal differences in surface 
heating doiiilliosnli . Therefore, we see in the Earth’s tropics 
on a small scale what tida lly locked terrestrial plane ts ex¬ 
perience on a global scale dShowman Polvanill201lL . The 
embedded reverse circulation cells, once formed, are present 
even for faster rotations IFig. 1171) . 

The abs e nce o f a state 1 circulation cell in the model of 
lEdson et aid d201ll ') and the presence of a second extended 
circulation cell per hemisphere already for Prot = 10 days 
can be explained by the dominance of the extra tropical 
Rossby waves in their model. On Earth, the second pair of 
circulation cells, the Ferrell cells, are located polewards of 
the direct tropical circulation cells and are a result of eddy 
circulation and are flanked by extra tropic al Rossby wav es 
(the subtropical and tropical jet, see, e.g. iHoltonI dl992ll f . 
Thus, we speculate that the second pair of circulation cells 
per hemisph e re for the planets with Prot ^ 10 days shown in 
lEdson et al.l (l201ll ') are likewise driven by eddy circulation 
and are associated with strong extra tropical Rossby waves 
that are stronger in their model than in our nominal model. 
The relative weakness of extra tropical Rossy waves in our 
model may also explain why the secondary circulation cells 


{state 2), polewards of the direct equatorial circulation cell, 
are much weaker in the Prot = 10 days simulation of an 
Earth-size pla net (Fig. 1171) than for the same set of planet 
parameters in lEdson et al.l il201ll l (their Eig. 4i). 

While for the transition from state 0 to state 1, a link 
can be found between circulation and standing (tropical) 
Rossby waves, no such link can be established at first glance 
for the other circulation state transitions with rotation pe¬ 
riod when inspecting Table [T] This is surprising because the 
Lp/Rp ~ 1 and Xp/Rp « 0.5 transition regions were previ¬ 
ously identified with suppression of the uprising direct circu¬ 
lation branch fSection 15.2.21) . Indeed, the equatorial direct 
circulation cell is very weak and also vertically suppressed in 
state 3. This state, however, is only established for relatively 
fast rotation (Prot ^ 4 days). 

A closer inspection of Figure 1181 top-left panel, shows 
for the Prot = 6 days simulation is illuminating: Although 
the direct equatorial circulation cell is not fully suppressed, 
it appears to be fragmented into a lower segment close to 
the surface between p — 1000 — 800 mbar and an upper seg¬ 
ment at p = 500 — 100 mbar by the embedded circulation 
cell. The disappearance of upwelling signatures at the top 
of the atmosphere over the substellar point in Section 15. 2.2 [ 
thus, would not indicate the complete disappearance of di¬ 
rect circulation but instead its fragmentation. As a conse- 
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Figure 18. Meridional mass flux stream function ^ in units of 10® kg/s for Rp = ^Rparth showing the four possible circulation states 
and the general reduction of circulation strength with faster rotation. Note that positive values indicate clockwise and negative values 
counterclockwise circulation. Contour intervals are 25 X 10® kg/s. Arrows indicate the upper segment of the direct circulation cells. 


quence, only a part of rising air masses over the substellar 
point is carried high aloft, whereas the other part is trans¬ 
ported close to the surface. Condensation of volatile species 
at high altitudes after evaporation at the surface, and thus 
cloud formation, may be suppressed in such a circulation 
state. We find the fragmentation to be the strongest for the 
Rp = 2REarth planet (Figure [TS]). It makes the identihca- 
tion of distinct circulation cells very difficult. It is doubtful if 
a closed transport cycle of condensables, like the water cycle, 
can be established in such fragmented circulation states. 

At the very least, the onset of circulation state 2 with a 
secondary full circulation cell per hemisphere should depend 
on rotation period and planet size, if one naively assumes 
that secondary circulation cells form once the direct circula¬ 
tion cells have shrunk sufficie ntly in latitudinal extent. Ac¬ 
cording to lHeld &: Houl (ll98Cll l. the latitudinal extent of the 
direct cell is 


V sn^Rl ) 


(25) 


Instead, we find that the transition to state 2 is shifted 
to smaller rotation periods with increasing planet size (Ta¬ 
ble 0. Equation!^ also predicts that the Hadley cell extent 
vh is smaller than 90° already for Prot = 100 days and for 
any planet size. According to these results, we should not see 
a state 0 circulation state in our parameter study. However, 
Equation [2S] is only truly applicable to direct unperturbed 
circulation cells, only for small Coriolis forces and was devel¬ 
oped to describe the circulation regime of the Earth, where 
the temperature gradient is between equator and pole and 
not between day side and night side. It is, thus, not really 
surprising that this simple relation fails to capture the evo¬ 
lution of circulation cells in our study. 

Direct circulation cells not only lose in strength with 
fast planet rotation and become fragmented, the cells can be¬ 
come almost fully suppressed in circulation state 3 IFig. 1181 
lowest panels, and 1 171 top-left panel): They appear ’pushed’ 
down towards the surface by the equatorial superrotation. 


where A6p-eq is the fractional change in potential tempera¬ 
ture from pole to equator. This relation predicts that the ex¬ 
tent of the direct Hadley circulation cell shrinks with planet 
size and angular velocity. Thus, one would expect a transi¬ 
tion to state 2 with more circulation cells to appear at slower 
rotations for a large planet when compared to a small planet. 


The suppression of the direct circulation starts at 
Prot ss 4 days for all planet sizes (Fig. 1161 left panel), 
whereas secondary circulation cells gather strength and even 
start to dominate (Fig. 1161 right panel) for large terres¬ 
trial planets {Rp = 1.25 — 1.75REarth)- See also Figure [TSl 
bottom panels, where we have added arrows for the better 
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identification of the direct circulation cell remnantfl Inter¬ 
estingly, the transition to circulation state 3 not only sup¬ 
presses direct circulation, it also stabilizes the other circu¬ 
lation cells against fragmentation for Rp — 2REarth (Fig- 
ure m bottom panels). A stable circulation state 3 may 
make the stable transport cycle of volatiles again possible. 
The transport, however, has to occur - in the absence of 
direct circulation - via the secondary and tertiary cells. 

Apparently, full 3D planet atmosphere models are re¬ 
quired to track the evolution of circulati on, which is a lready 
known from the Earth climate regime (lHoltonlll992li . But 
even then we need to be wary of differences in the climate 
dynamics - like the relative strength of the tropical and extra 
tropical Rossby wave - that may lead to drastic deviations 
between the assumed circulation state. 

The suppression of direct circulation for fast rotating 
terrestrial planets that we rep ort here is indeed ve ry different 
from the results obtained bv lEdson et al.l (l201llf . where di¬ 
rect circulation is always dominating. Even for Prot = 1 day 
with a mixed state between high latitude and equatorial 
westerly jets, the direct circulation cells are not suppressed 
but flanked by two jets at ± 10°. The wind speeds are, 
however, much slower (u « 100 m/s) than the ones found in 
our simulations (u ~ 300 m/s). iMerlis fc Schneider! (l2010l j 
don’t have superrotation but high latitude westerly in their 
Prot = 1 day simulation. Their direct equatorial and sec¬ 
ondary circulation cells are of similar strength. All three 
studies are in different climate states for an Earth-sized 
planet with Prot = 1 days, which explains why they dis¬ 
agree in the details of the circulation state. At least they 
agree that a fast rotating tidally locked Earth-size planet 
has three circulation cells per hemisphere. 


7 TEMPERATURES AND HABITABILITY 

Therefore, we have shown once again the merit of our de¬ 
tailed and coherent parameter space investigation with re¬ 
gards to standing Rossby waves. We can tie our improved 
understanding of the Rossby wave climate regimes derived 
in the previous section to changes in the circulation states. 
We can even coherently explain differences in climate states 
between different models. 

Changes in circulation and in climate dynamics due to 
different standing planetary wave configurations have an in¬ 
fluence on surface temperatures and thus, potentially, on 
habitability. Thus, we discuss in the following the evolution 
of surface temperatures with planet size and rotation period 
in the light of the previously gained insights. 

7.1 Surface temperature dependency on planet 
size 

First of all, bigger planets tend to have larger surface tem¬ 
perature gradients in our study: they have a hotter substellar 
point and cooler night sides than smaller planets IFigure fT^ . 

^ As a reminder, direct circulation is counter-clockwise in the 
South hemisphere (negative latitudes) and clockwise in the North 
hemisphere (positive latitudes), for the secondary circulation the 
circulation direction is reversed (see Figure [TTll . 


Furthermore, the average surface temperature is higher for 
small planets compared to big planets. These results are 
surprising as the circulation strength appears to be stronger 
for larger planets (Figure [T6l left panel) and also the wind 
speeds are higher (Figure [S]). Both effects would suggest a 
better temperature mixing and thus a smaller temperature 
gradient for large planets. However, the surface temperature 
evolution not only depends on wind speeds and circulation, 
that is, on planet atmosphere dynamics, but also on the 
stellar irradiation. In the context of stellar irradiation, the 
eternal night side in our model always cools down and the 
eternal day side is always irradiated. The efficiency of the ra¬ 
diative cooling/heating with planet size can be estimated via 
Equation [15] The efhciency increases with increasing planet 
size and thus increasing surface gravity g. The dynamical 
time scale Tdyn = U/Rp (Section [3|, on the other hand, is 
expected to also increase with Rp, at least for slow rota¬ 
tion, because we would expect U ~ Rip. Thus, the increase 
in radiation efficiency is naively expected to be balanced 
by increased dynamical efhciency. The comparison between 
zonal wind maxima in Section [5. 1.1 1 indicates that the wind 
speeds indeed increase as expected with planet size. The re¬ 
lations discussed in Section FS.l.ll however, are only valid for 
zonal wind velocities above the planetary boundary layer, 
which is true for the discussed maximum zonal wind speeds. 
For surface temperatures, however, surface friction can no 
longer be neglected. Furthermore, it already has been shown 
in Section 15.1.11 that U ~ Rip is only a valid assumption 
for slow rotation with negligible Coriolis force. Therefore, in 
the fast rotation regime, the increase in wind velocities with 
planet size is smaller than estimated by the scale analysis. 
For these two reasons, the dynamical efhciency of surface 
how does increase with planet size but not as fast as the 
efhciency of radiative forcing. Consequently, larger planets 
have cooler night sides and hotter day sides. 

7.2 Surface temperature dependency on planet 
rotation period 

The dependency of surface temperatures with rotation is 
more complex but follows tendencies that have been iden- 
tihed earlier with respect to horizontal how (Section 15.21) 
and circulation (Section O. Slow rotating planets (Prot = 
30 — 100 days) are well mixed by the direct equatorial cell 
{state 0) and, thus, show the smallest surface temperature 
gradient (Figure fT^ upper panels). With faster rotation and 
thus diminishing direct circulation strength (Figure fTOl left 
panel), the night side temperatures drop more rapidly than 
the substellar point temperatures, which leads generally to 
lower average surface temperatures (Figure [T9l lower panel). 
Interestingly, the day side surface temperatures also decrease 
slightly with decreasing circulation strength. This suggests 
that the night side branch of the circulation that heats the 
night side via transport from the hot substellar point, is 
more strongly diminished than the day side branch that 
cools the substellar point via upwelling. The difference can 
be explained by taking into account surface friction that 
diminishes the surface flow and thus acts primarily on the 
lower branch of the circulation with fl ow from the n i ght si de 
towards the substellar point (see also lCarone et aP 1120141) ). 
The longer air masses remain on the night side, the more 
heat can be radiated away, the cooler the air becomes. The 
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Figure 19. Minimum surface temperatures (top-left panel), maximum surface temperatures (top-right panel) and average surface 
temperatures (lower panel) versus rotation period and planet size for Earth-like irradiation and an Earth-like atmosphere. The dashed 
gray line marks T = 273 K, the solid gray line marks T = 373 K. 


vertical uprising branch directly over the substellar point is, 
on the other hand, only indirectly affected by surface fric¬ 
tion. The net effect is an overall cooling of the surface. 

The transition from a state with only direct circulation 
{state 0) to a state with direct circulation and embedded re¬ 
verse circulation cells {state 1) can be identified in Figure [T^ 
by a change in the cooling rate of average surface tempera¬ 
tures between Prot = 20 — 34 days, depending on planet size, 
respectively (Compare with Table [7|). The circulation tran¬ 
sition results in a strong decrease in minimum night side 
temperatures due to the reduced efficiency of heat trans¬ 
port from the day to the night side (Figure [TG] left panel). 
Interestingly, the substellar point surface temperature also 
cools down although this effect is less pronounced for larger 
planets. 

With even faster rotations, the investigated climate sim¬ 
ulations approach the Lr « 1 and Xr « 0.5 transition re¬ 
gion (Section 15.2.21) . We associated this region with frag¬ 
mentation of the uprising branch of the direct circulation 
over the substellar point and thus reduced upwelling over 
the substellar point. At the same time, secondary circula¬ 
tion is still very weak. Thus, it is puzzling that we find in 
this region apparently a warmer night side and average sur¬ 


face temperatures {Prot = 8—12 days for Rp = IRsarth 
and Prot = 9 — 16 days for i?p = 2REarth, compare Table [5] 
with Figure [Toll . The fragmentation of the direct circulation 
cells into a lower and an upper segment may explain the 
warmer night sides. We already speculated in Section [6] that 
air masses are transported closer to the surface in the lower 
segment, which would represent a ’short cut’ compared to a 
circulation state with a fully vertically extended circulation 
cell. The lower circulation segment may thus heat the night 
side surface efficiently. 

It is puzzling though that day side temperatures re¬ 
main constant for large planets {Rp > l-75REarth), until 
the Lr/Rp « 0.5 transition with Prot ~ 4 days is reached 
(Table[5}. For small planets {Rp ^ lAHREarth), the day side 
temperatures cool down with faster rotation, as expected, 
because direct circulation strength decreases with faster ro¬ 
tation. The days side cooling stops at the Lr/Rp « 0.5 
transition region {Prot = 5 — 6), where the simulations reach 
another atmosphere dynamics regime. 

The rotation regime with Lr/Rp ^ 0.5 brings a rise 
in surface temperatures for all planet sizes, where the sub¬ 
stellar point temperatures become very hot and the night 
side is moderately warmed. The strong increase in day side 
surface temperatures can be explained by the suppression of 
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Figure 20. Zonally and temporally (1000 days) averaged potential temperature (upper panels) and temperature in K for an Earth and 
2 Earth size planet and rotation period Prot = 1 days. Contour intervals are 10 K. Arrows indicate location of thermal inversions. 


the equatorial direct circulation cell and thus upwelling that 
would normally cool down the substellar point efficiently. At 
the same time, secondary circulation cells gain in strength 
and manage to transport heat from higher latitudes at the 
day side towards the night side - albeit not as efficient as the 
direct circulation in the slow rotation regime Prot is 30 days. 
This explains the moderate warming of the night side sur¬ 
face. 


The fragmented circulation for Prot = 2 — 6 days on 
the 2 Earth radii planet is reflected in fluctuating night side 
minimum temperatures (Figure fTSll . Apparently, the night 
side is even more strongly affected by this strange circulation 
than the day side. 


In summary, only the circulation transition from state 
0 to state 1 and the fragmented circulation of the 2 
Earth-radii planet correspond to changes in surface tem¬ 
peratures. Fragmented circulation also mainly affects the 
night side. Climate dynamics transitions brought about by 
changes in standing Rossby waves, on the other hand, play 
again a strong role in shaping surface temperatures via 
ch anges in the horizon tal flow. We have already discussed 
in ICarone et al.l (|2014| ) that our surface temperatures are 
at the night sides cons iderably warmer th an i n the comple x 
Earth model studies of I Edson et_alj j[20n ) an OS hl ll2003D . 
while they agree roughly with Joshi et al (IiOOtI ). The latter 


also used a model with simplified forcing. We will show in 
an upcoming follow-up study, that it i s possible to c hange 
the temperature forcing prescription of ICarone et al.l (l2014l ) 
to include more efficient night side cooling. We will discuss 
there how and if more efficient night side cooling changes 
the results derived in this study. 


7.3 Vertical temperature profile 

The surface temperatures are only a small part of the to¬ 
tal temperature structure in a planet’s atmosphere. In the 
following, we will thus discuss the complete vertical temper- 
atu re structure. _ 

lEdson et al.l il201ll ) showed that the slope of the isen- 
tropes in potential temperature between the tropics and the 
mid-latitudes is much gentler for a tidally locked Earth-size 
planet with Prot = 1 days than in the dry Earth control run. 
I Edson et al.l (l201ll ) attributed that to a standing Rossby 
wave that we have identified in this study as a standing ex¬ 
tra tropical Rossby wave. 

The vertical potential temperature profile for a tidally 
locked Earth-size planet with Prot = 1 day in our study (Fig- 
ur el20l) exhibits muc h steeper isentropes than in Figure 4b 
of lEdson et al.l (l201ll ). This difference is again evidence that 
our simulations are in a ver y different cl i mate dynamics 
regime than the simulations of lEdson et al.l (l201ll ): one that 
is dominated by standing tropical instead of extra tropical 
Rossby waves. Our results show, furthermore, that the isen¬ 
tropes are steeper for a larger planet {Rp — 2REarth)- These 
results are expected as the steepness of the isentropes are 
also determined by the strength of the Coriolis force, which 
is stronger on a larger planet. Furthermore, we see evidence 
of temperature inversion at the top of the atmosphere for 
this very fast rotation, that is, a local increase in tempera¬ 
ture with height instead of a decrease. 

We will demonstrate, in the following, that two different 
types of upper atmosphere thermal inversion can be found 
in our simulation results. One type of thermal inversion de¬ 
velops because of chevrons around the very fast superrotat- 
ing equatorial jet. The other type was already mentioned 
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Figure 22. Zonally and temporally (1000 days) averaged temperature in K for a 2 Earth size planet around the Ln/Rp = 1 and 
Xp/Rp = 0.5 region, that is, for Prot = 16,13,9,7,4, and 3 days, respectively. Contour intervals are 2 K. Arrows indicate location of 
thermal inversions. 


in Sections 15.2.11 and 15.2.21 Thermal inversions form at the 
location of the previously identihed local temperature max¬ 
ima over the cyclonic vortices at p = 225 mbar, which are 
the result of adiabatic heating of air masses flowing down¬ 
wards towards region of reduced geopotential height. We 
found in Section |6l however, that thermal inversion should 
only be present for Rossby climate state regimes between 
Xp/Rp < 0.5 and A ^ 1. 

When inspecting the horizontal cross section at the 
top of the atmosphere for very fast rotating planets with 
Prot = 1 days (Figure [T4ll . it becomes indeed clear that the 
temperature inversion can be linked to the ’chevron’ that 
the superrotating equatorial jet develops, in particular, on 
larger planets. The ’chevron’ for the fast rotating large plan¬ 
ets links our study to atmosphere mo dels of hot Jupi t ers th at 
also show this ’chevron’ shape (e.g. iKataria et al.l ll2015l) ). 
The upper atmosphere temperature inversions associated 
with the chevron does not change in size and altitude with 
larger planet size but expand horizontally towards higher 
latitudes. 


We hnd correspondingly the second type of temperature 
inversions for slower rotations (Figure [2l1l . We could also 
indeed verify that they start to develop at the Xp/Rp = 1 
transition region and are not present for very slow rotation 
with Prot >>45 days (not shown). These type of tempera¬ 
ture inversions drop in altitude and increase in vertical ex¬ 
tent for larger planet sizes, because of the steeper isentrope 
slopes in the mid-atmosphere {p = 525 mbar) below the 
cyclonic vortices, in which the upper atmosphere flow gets 
adiabatically heated. The upper atmosphere flow can ’fall’ 
deeper and farther towards the poles. 

Furthermore, we expect these temperature inversions to 
vanish, when the rotation period approaches the Lp/Rp =1 
and Xp/Rp = 0.5 transition region (Figure [22ll . at the same 
time as the adiabatic heating over the vortices diminishes 
and disappears (Section l5.2.2ll . 

Indeed, we can show that the temperature inversions 
reduce in vertical extent and are shifted towards the very top 
of the atmosphere for faster rotation, despite the steepening 
of the isentropes that should favour adiabatic heating and 
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Figure 21. Zonally and temporally (1000 days) averaged tem¬ 
perature in K for an Earth (upper panel) and 2 Earth size planet 
(lower panel). In both cases: Prot = 16 days. Contour intervals 
are 2 K. Black arrows mark the location of the upper atmosphere 
thermal inversion. 


thus increase adiabatic warming of the upper air flow. We 
focus again on simulations for a large planet {Rp = ^Rparth) 
with Prot = 3 — 16 days (Figure I22|l . The shift towards 
the sponge layer explains why we saw the thermal maxima 
disappearing already at Prot = 9 days, when we inspected 
the horizontal temperature map at p = 225 mbar (Figure^). 

The weakening of the temperature inversion can again 
be explained by suppressed upwelling due to circulation 
cell fragmentation (Figure [TS] top and middle panels). For 
Prot = 4 day, the feature disappears completely. With even 
faster rotation, the chevron shape develops and forms a new 
temperature inversions. It appears first at the equator and 
extends towards higher latitudes with faster rotation (Fig¬ 
ure lower panels and Figure O. We have to caution 
that the equatorial jet, the chevron and thus the associated 
thermal inversion are very close to the sponger layer. Fur¬ 
thermore, we don’t resolve the upper atmosphere very well 
due to our strong focus on the investigation of dynamics in 
the troposphere. It will have to be investigated, if this type 
of thermal inversion extends to higher altitudes in a model 
with a better upper atmosphere resolution. 


8 OBSERVATIONAL DISCRIMINATION 

BETWEEN DIFFERENT CLIMATE STATES 

It has been shown that for Lp/Rp ^ 1 and Xr/Rp ^ 0.5 dif¬ 
ferent climate states are possible for the same set of param¬ 
eters. These degeneracies in climate states should be taken 
into account at least for tidally locked terrestrial planets at 
the inner edge of the habitable zone around M dwarfs (Sec¬ 
tion [2]4]). We will now highlight several simulation results 
that provide insights into the observational discrimination 
between different climate states. 

In a climate state dominated by a standing tropical 
Rossby wave, the hottest point is displaced eastward with 


respect to the substellar point by strong equatorial superro - 
tation as already predicted bv IShowman fc Guillotl 1120021) . 
We And for the Prot = 1 day simulations that the hottest 
point is indeed shifted eastward - as expected (Figure [14]). 
The eastward displacement of the hot atmosphere region 
was inferred for the transiting hot Jupiter HD 189733 b from 
Spitzer lighcurves. The orbital phase variations showed that 
the maximum of the planet’s phase curve and thus the max¬ 
imum in planetar y thermal emission o c curre d prior to the 
secondary eclipse dKnutson et ahllioOTl . 1200^ . 

The hottest point is, however, not shifted eastward but 
westward with respect to the substellar point in climate 
states with a strong extra tropical Rossby wave component 
(Figure nn top and middle panel and Figure [TS] left panel). 
Such a climate state would thus show up in the planet’s 
orbital phase curve because the maximum would occur after 
the secondary eclipse. 

The Prot = 3 days simulation for Rp = 1.25REarth 
shows another mixed climate state, one in which standing 
tropical and extra tropical Rossby waves contribute in ap¬ 
proximately equal measures (Figure fT^ left). The result is a 
hot region eastward and westward from the substellar point 
(Figure 1131 l eft). The temperatu re m ap is reminiscent o f 
the results of iKataria et al.f (l2014l ) and [Lewis et al.l ll2010l) . 
who also find mixed states with high latitude jets and a rela¬ 
tive strong equatorial superrotation in simulations with solar 
metallicity for the Super-Earth GJ 1214b and hot Neptune 
GJ 436b. They assumed, however, gas giant composition, 
that is, the absence of a solid surface. A climate state with 
a smeared out hot atmosphere region around the substellar 
point would show up as a broadeni ng of the planets the rmal 
emission around secondary eclipse (iKataria et al.l[2015l ). 

The relative low atmosphere temperature (T ~ 300 K) 
and small horizontal temperature differences {ATmax ~ 
50 K) discussed in this study, preclude currently the observa¬ 
tional identification of different climate states for a transit¬ 
ing habitable planet. It may be possible, however, to discrim¬ 
inate between different climate states for hot Super-Earths. 
Future work will, therefore, extend this study to scenarios 
with stronger thermal forcing. 

Interestingly, the Kepler orbital phase curve of the hot 
Jupiter Kepler 7b shows indeed westward displacement of 
the maxi mum in planetary flu x with respect to the substel¬ 
lar point (jPemorv et al.l[2013l ). Thermal emission was ruled 
out by the authors because superrotation and thus an east¬ 
ward displacement of the hottest atmosphere was expected. 
The westward displacement was instead attributed to re¬ 
flected stellar light from inhomogeneous planet cloud cov¬ 
erage. It is too early to infer that Kepler 7b may be in a 
climate state with a strong extra tropical Rossby wave com¬ 
ponent. It would also be a puzzling result, because extra 
tropical Rossby wave s are only expected for sm all planets 
with a solid surface ll Showman fc Polvanil I 2 OI 1 I ) . We find, 
indeed, the strongest signals of extra tropical Rossby wave 
components for the smallest planets in our study (Figures 

13 and M with Rp — ^Rparth and \.25REarth^' 

On the other hand, iKataria et al.l |201^ showed that 
hot Jupiter climate models that impose large frictional drag 
exhibit an equatorial superrotating jet and two additional 
westerly jets at mid-latitude - just like in the climate states 
with a mixture of tropical and extra tropical Rossby waves 
identified in this study. Thus, large Super Earths with a 
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solid surface boundary dZalucha et aljlioisf) and gas plan- 
ets that experience frictional drag inside their atmospheres 
l|R,auscher fc Menoul l2013l : iKataria et alJ I201 !tI 'I may show 
similar climate state degeneracies. In an upcoming study, 
we will show that the appearance of a strong extra tropical 
Rossby wave component in large scale circulation is strongly 
tied to frictional drag at the surface boundary. 


9 SUMMARY 

We found that standing tropical and extra tropical Rossby 
waves are responsible for different climate regime phases on 
tidally locked habitable planets in this very thorough pa¬ 
rameter space investigation that covers Prot = 1 — 100 days 
and Rp = 1 — 2REarth- Thus, we could tie climate state 
phase transitions to transitions in the meridional Rossby 
wave numbers Xp/Rp and Lp/Rp for the tropical and ex- 
tratropical Rossby wave number, respectively. 

• The formation of a standing tropical Rossby 
wave at Xp/Rp = 1 brings about fast er zonal wind speeds, 
which is in line with the explanation of IShowman fc Polvanil 
ll201lli that the coupling between standing tropical Rossby 
waves and Kelvin waves transport angular momentum from 
high to lower latitudes. Furthermore, cyclonic vortices gain 
in strength and zonal wind maxima at mid-latitudes appear 
that get more and more confined towards the equator with 
increasing rotation. The vortices are cold features in the 
lower and middle atmosphere and thus have a low geopo¬ 
tential height. If upwelling flow over the hot substellar point 
with large geopotential height ’falls into’ these vortices, the 
flow is adiabatically heated which leads to the formation 
of temperature inversions in the vertical temperature struc¬ 
ture. Furthermore, small embedded reverse circulation cells 
form inside the direct circulation cells. The change in circu¬ 
lation state leads to a decrease in night side surface temper¬ 
atures. 

• The Lp/Rp = 1 transition region and Xp/Rp = 0.5 
transition region almost overlap in rotation period. Because 
our model favours a climate state dominated by tropical 
Rossby waves, the Xp/Rp = 0.5 transition is vastly 
more important in our model. Once the transition is 
crossed, zonal wind speed - generally - increases strongly 
with faster rotation. The transition in tropical Rossby wave 
number is also associated with beginning fragmentation of 
the direct circulation cell into a lower and upper segment, 
which leads to a decrease in upwelling at the top of the at¬ 
mosphere. As one consequence, the upper atmosphere tem¬ 
perature inversions diminish. At the same time, the lower di¬ 
rect circulation segment allows to warm the night side more 
efficiently, which leads to an increase in night side surface 
temperatures. 

• In contrast to our simulation results, the Lp/Rp = 1 
trans ition is more important in the model of lEdson et al.l 
1 I 2 OIII) . The authors report abrupt transition to lower zonal 
wind speeds. They, furthermore, describe the transition from 
a climate state with a broad equatorial zonal jet to a state 
with two zonal j ets at higher l atitud es. Thus, we conclude 
that the model of lEdson et al.l (1201 ill assumes a second pos¬ 
sible climate state in this rotation period region dominated 
by the extra tropical Rossby wave. 


• With transition to Lp/Rp ^ 0.5, both the tropical 
and extra tropical Rossby wave can ’fit’ on the planet. And 
indeed, we find several examples of climate states that show 
clearly the simultaneous presence of the tropical and extra 
tropical Rossby waves. Even for Prot = 1 days, there still 
exist remnants of a standing extra tropical Rossby wave 
that leads to broadening of the equatorial jet. In general, 
however, the tropical Rossby wave dominates even in this 
climate state regime - at least in our model. 

• The dominance of tropical Rossby waves explains the 
formation of a single strong equatorial superrotating jet with 
very fast wind speeds of up to 300 m/s with very fast ro¬ 
tations Prot ^ 1.5 days. We have to caution, however, that 
the exact value of the wind speed should be taken with a 
grain of salt. The superrotating jet are pushed with faster 
rotations into our sponge layer, whose friction efficiency 
needed to be increased to ensure numerical stability. The 
physical mechanism behind the formation of fast superrota¬ 
tion is, however, already firmly esta blished from the atm o- 
sphere modelling of hot Jupiters fe.g. IKataria et ^ ll2015h '). 
Furthermore, we can explicitly demonstrate the meridional 
tilt of Rossby wave gyres arisi ng from the standing t ropica l 
Rossby wave as proposed by IShowman fc Polvanil ll201llf . 
Thus, we argue that while we are uncertain about the exact 
value of the wind speeds, the qualitative assessment about 
the formation of an equatorial superrotating jet with Rossby 
wave number Xp/Rp ^ 0.5 is sound. 

• The fast equatorial superrotation for Prot ^ 1-5 days 
suppresses the direct equatorial circulation cell practically 
completely. Circulation is instead provided by the secondary 
circulation cells whose strength increases with faster rota¬ 
tion. The suppression of the direct circulation is in stark 
contra st to Earth circulation studies dNavarra fc Boccalettil 
1 I 2 OO 2 IB and that of lEdson et al.l (1201 il l , where the direct cir¬ 
culation cell is always the strongest cell. These studies don’t 
show, however, strong equatorial superrotation, which ap¬ 
pears to vertically push the already fragmented direct circu¬ 
lation cell by the embedded reverse circulation cells towards 
the surface. Apparently, superrotation and direct circulation 
don’t mix very well. The suppression of the direct circulation 
prevents the substellar point at the surface to cool down via 
upwelling: Thus, the substellar point surface temperatures 
reach high temperatures. 

Furthermore, it was found that larger planets show in 
general larger wind speeds, steeper isentropes, more pro¬ 
nounced warm upper atmosphere temperature inversions, 
warmer substellar point and cooler night side surface tem¬ 
peratures compared to an Earth-size planet with the same 
rotation period. Furthermore, the largest planet {Rp — 
2REarth) exhibits strongly fragmented circulation between 
Prot = 2 — 6 days in contrast to smaller planets. Appar¬ 
ently, circulation and thus potentially cloud formation may 
be disrupted in this rotation regime. 

In addition, large very fast rotating planets {Rp = 
2REarth and Prot = 1 day) develop features that are al¬ 
ready known from the atmosphere modelling of hot Jup iters, 
Neptunes and Super -Earths (e.g. IKataria et al.l (l2015ll and 
IZalucha et al.l (l2013h f: Counter rotating jets at higher lati¬ 
tudes that flank the very fast equatorial superrotating jets, 
and warm chevrons at the substellar point. Apparently, large 
planet size and fast rotation are already sufficient for the for- 
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mation of such features, even though we keep at all times 
temperate Earth-like forcing. 


10 CONCLUSIONS & OUTLOOK 

We performed a parameter study with unprecedented thor¬ 
oughness in planet rotation period {Prot = 1 — 100 days) 
and planet size {Rp — l — 2REarth) to coherently investigate 
climate state changes of tidally locked habitable terrestrial 
planets around M dwarfs. We were not only ca pable of iden¬ 
tifying mor e climate state trans i tion s tates as lEdson et al.l 
(12011 1 and iMerlis fc Schneider! (l201Cll l. we also confirmed 
the existence of multiple climate states in the fast rotation 
regime {Prot < 12 days for Rp = 2REarth)- The existence 
of multiple c l imate state solutions was already indicated by 
lEdson et al.l ll201lll . who speculated about the involvement 
of the extra tropical Rossby wave. We were, however, capa¬ 
ble of identifying the possible formation of either the stand¬ 
ing tropical or the extra tropical Rossby waves as the origin 
of the multiple climate states. Our model apparently favours 
the first state w hich exhibits fast e quatorial superrotation 
and the model of lEdson et al.l ll201lli favours apparently the 
second state with two high to mid-latitude westerly jets with 
slower wind speeds. 

We attribute the differences be t ween the results from 
our model and that of lEdson et al.l (I2OIII ') to the relevant 
strength of the tropic al Rossby wave in our model in con¬ 
trast to the model of lEdson et al.l (I2OIIII . where the extra 
tropical Rossby wave dominates. We will show in an up¬ 
coming paper that an increase of surface friction efficiency 
in our model can also produce climates that are dominate d 
by extra tropical Rossby waves like in lEdson et al. I (l201lli . 
Reduction of surface friction time scales, on the other hand, 
increases the strength of tropical Rossby waves and leads 
to even faster equatorial superrotation. The coherent dis¬ 
cussion of climate state regimes and Rossby wave properties 
allowed us, thus, to identify sources of discrepancies between 
different models. 

We found, in addition, another region of climate state 
degeneracies that lEdson et al.l (I 2 OIIII missed: For very fast 
rotations {Prot ^ 12 days for Rp = 2REarth), a third possi¬ 
ble climate state solution appears, because the tropical and 
extra tropical Rossby wave can both form standing waves 
and allow for the existence of a mixed climate state with 
consequently three zonal wind jets. 

Interestingly, the upper atmosphere’s hot spot is east¬ 
ward shifted with respect to the substellar point in the first 
state, westward shifted in the second state and the third 
state shows a longitudinal ’smearing’ of the spot across the 
substellar point. Such climate states could potentially be 
observed in the orbital phase curve of transiting planet, be¬ 
cause the maximum would occur prior, after or ’smeared 
around’ the secondary eclipse. The relative low atmosphere 
temperature (T « 300 K) and small horizontal temperature 
differences {ATmax ~ 50 K), however, make the observa¬ 
tional identification of different climate states for a transit¬ 
ing habitable planet currently infeasible. It may be possible, 
however, to discriminate between different climate states for 
hot Super-Earths from the planet orbital phase shift with re¬ 
spect to the secondary eclipse. Future work will, therefore, 
extend this study to scenarios with stronger thermal forcing. 


How much are the investigated climate states of rel¬ 
evance for habitable planets? The planet sizes between 
IRsarth and 2REarth are indeed relevant as this planet mass 
and size range is compatible with rocky terrestrial Super- 
Earth planets dMadhusudhan et al.ll2012 ). The relevance of 
atmosphere dynamics for short orbital periods Prot ^ 8 days 
is more difficult to assess. The convent i onal definition of 
the ha bitable zone dKasting et al.l dl993lllKoDparaDU et al.l 
d2013tl ') puts the inner edge of the habitable zone for the 
smallest M dwarf star (M« = 0.2Msun) at 0.08 AU which 
corresponds to a rotation period of 18 days for the fastest 
rotating habitable planet. TVang et aP d2014l i showed that 
the habitable zone can be expanded towards the star due 
to cloud formation at th e substellar point , which increases 
the albedo. Furthermore. IZsom et al] d2013l ') showed that the 
habitable zone can be pushed even closer to the host star, 
if the terrestrial planet is assumed to be relatively dry. The 
inner edge is located at 0.04 AU for M, = 0.2Msun, which 
corresponds to Prot =6.5 days. With this new definition of 
the inner habitable zone, our results become very relevant: 
A large Rp = 2REarth planet may, for example, exhibit 
fragmented circulation cells for Prot = 6 days, which raises 
immediately the question, if water clouds can form. 

The very fast rotation regime Prot ^ 1-5 days provides, 
in any case, already a link towards the study of atmosphere 
dynamics on gas planets: We have shown that the largest 
planet in our study develops fast equatorial superrotation 
and a temperature chevron that are atmosphere dynamics 
features well known from hot Jupiters. At the same time, 
even our largest planet still has remnants of extra tropical 
Rossby waves that are la cking on hot Jupiters due to their 
lack of a solid boundary JShowman fc Polvanill2011fl . Thus, 
our study provides already at this stage insights about the 
transition of climate atmosphere dynamics from small hab¬ 
itable terrestrial planets to larger Mini-Neptunes. This tran¬ 
sition region is of great interest for Super-Earth planets that 
sho w no clear-cut distinct i on be t ween rocky and g a s plan ets 
('e.g. iMadhusudhan et al.l (l2012l i lRogers fc Seager (l2010l 'li. 

As already mentioned, we will strengthen in upcom¬ 
ing follow-up studies the connection between extra tropical 
Rossby waves and the presence of a solid boundary. The 
westward displacement of the upper atmosphere’s hot spot 
with respect to the substellar point may then reveal the pres¬ 
ence of a solid surface or other sources of frictional drag. 
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